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INTRODUCTION 

The  purpose  of  this  document  is  to  provide  an  abbreviated  user's  guide  for 
the  current  version  of  the  Rotocraft  Flight  Simulation  Program  C81  as  installed 
at  NAVAIRDEVCEN .  This  document  is  to  be  used  as  a  quick  guide  to  the  program 
capabilities  of  C81  in  the  area  of  Flight  Dynamics  analysis.  It  also  describes 
the  problems  encountered  with  the  current  version  of  the  program,  designated 
AGAP80,  as  delivered  to  this  Center.  It  should  be  noted  that  the  problems  that 
are  described  herein  are  associated  with  the  copy  of  the  CDC  converted  program 
that  is  installed  on  the  CDC  CYBER  175/CDC  6600  Computer  System  at  NAVAIRDEVCEN 
and  are  not  considered  universal. 

The  Rotorcraft  Flight  Simulation  Program  C81  has  been  developed  and  under 
continual  update/revision  since  the  early  1960's  by  Bell  Helicopter  Textron  under 
contract  with  the  Applied  Technology  Laboratory,  U.S.  Army  Research  and 
Technology  Laboratories  (AVRADCOM) ,  Ft.  Eustis,  VA.  All  versions  of  the  program 
have  been  developed  for  usage  on  an  IBM  computer  system.  Early  in  1980  the  Naval 
Air  Development  Center  and  the  NASA  Langley  Research  Center  expressed  interest  in 
a  CDC  conversion  of  the  program  for  their  CDC  CYBER  computer  systems.  A 
successful  conversion  of  the  program  was  completed  in  June,  1980  by  Mr.  Robert 
Tomaine  of  the  Army  Structures  Laboratory  and  a  magnetic  tape  of  the  program 
source  was  delivered  to  NAVAIRDEVCEN  in  July,  1980. 

This  conversion  was  done  to  the  most  current  version  of  the  computer  program 
and  is  referred  to  as  AGAP80. 

After  many  installation  problems  were  overcome,  the  program  was  installed  on 
the  CDC  CYBER  175/CDC  6600  Computer  System  at  NAVAIRDEVCEN  in  March,  1981.  The 
period  from  March  through  September,  1981  was  devoted  to  validation  of  the  CDC 
conversion  on  NAVAIRDEVCEN 's  computer  system  and  understanding  of  the  program 
capabilities  and  applications. 

A  complete  input  data  deck  for  an  AH-1G  helicopter  was  provided  with  the 
source  tape  as  a  sample  case  for  validating  the  program  at  NAVAIRDEVCEN.  This 
AH-1G  data  deck  is  used  throughout  this  report  to  illustrate  the  program 
capabilities. 

An  effort  is  made  throughout  the  description  of  the  program  capabilities  to 
provide  the  user  with  a  description  of  the  inputs  required  to  perform  that 
particular  program  operation  and  the  output  data  that  resulted.  The  user  is 
referred  to  the  user's  manual  for  the  AGAP80  version  (reference  1)  for  the 
information  necessary  for  assembling  a  complete  input  data  deck  and  to 
successfully  execute  the  program. 

Two  program  operation  problems  were  identified  while  running  test  cases  to 
check  out  each  of  the  program  operational  capabilities.  These  two  were  the 
inability  to  perform  any  of  the  postprocessing  operations  of  GDAP80  and  the 
inability  to  compute  the  rotorcraft  characteristic  equation  roots  and  associated 
transfer  functions  in  the  rotorcraft  stability  analysis.  Both  of  these  problems 
are  associated  with  the  copy  of  the  AGAP80  version  that  is  installed  at 
NAVAIRDEVCEN  and  are  not  universal  problems  with  this  program.  The  cause  of  the 
postprocessing  problem  has  been  identified  and  is  relatively  easy  to  fix,  while 
the  eigenvalue  calculation  problem  is  still  under  investigation. 
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This  program  provides  the  capability  to  analytically  examine  the  flying 
qualities  characteristics  of  a  given  model  of  helicopter,  whether  it  be  single, 
tandem,  or  tilt  rotor.  The  input  data  can  be  minimal  for  a  rough  first-cut 
examination  of  a  preliminary  design,  or  extensive  for  a  simulation  of  an  existing 
aircraft  for  which  fuselage  aerodynamic,  airfoil,  aeroelastic  blade,  and  rotor 
induced  velocity  distribution  data  are  already  known  through  previous  testing. 

This  program  can  perform  a  trimmed  flight  analysis,  a  rotorcraft  stability 
analysis,  parameter  sweeps  of  trimmed  flight  conditions,  and  a  maneuver 
simulation.  It  also  can  retrieve  maneuver  time-history  data  stored  on  magnetic 
tape,  print  plots  of  time-history  data,  and  store  maneuver  time-history  data  on 
tape.  The  program  has  other  postprocessing  capabilities  which  are  performed  on 
time-history  data  which  would  be  of  interest  to  structural  analysis  engineers. 

All  of  these  operational  capabilities  are  discussed  in  this  report. 
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PROGRAM  CAPABILITIES 


AGAP80  OPERATIONS 

The  AGAP80  version  of  the  C81  program  is  capable  of  performing  the  following 
general  operations: 

1)  Compute  a  trimmed  flight  condition 

2)  Perform  a  rotorcraft  stability  analysis 

3)  Perform  parameter  sweeps  of  trim  conditions,  with  or  without  a 
rotorcraft  stability  analysis 

4)  Compute  a  maneuver  with  or  without  a  rotorcraft  stability  analysis 

5)  Retrieve  maneuver  time-history  data  stored  on  magnetic  tape 

These  five  operations  are  illustrated  in  the  flow  charts  given  in  figures  1 
through  8.  The  block  labelled  "OPERATIONS  ON  TIME-HISTORY  DATA"  represent  the 
execution  of  the  postprocessing  program,  GDAP80. 

Each  of  the  AGAP80  operations  or  combination  of  operations  is  controlled  by 
input  data.  Since  the  amount  of  data  to  be  read  depends  on  the  operation  or 
operations  desired,  a  data  deck  in  this  context  consists  of  a  message  card,  an 
"NPART"  card  telling  the  program  which  primary  operation  or  operations  to 
perform,  and  the  additional  data  necessary  to  perform  the  indicated  operations. 

In  some  cases  the  additional  data  are  contained  on  500  or  more  additional  cards, 
while  in  other  cases,  as  little  as  a  few  cards  of  additional  data  is  required  for 
the  AGAP80  input  deck. 

The  second  step  in  several  of  the  flow  charts  is  "CALCULATE  PROBLEM 
CONSTANTS."  In  each  operation  containing  this  step,  a  number  of  quantities  which 
remain  constant  throughout  the  performance  of  the  operations  must  be  defined 
using  the  input  data.  Performing  such  computations  drastically  reduces  the 
number  of  program  inputs  and  also  provides  program  flexibility  necessary  for 
incorporating  such  operations  as  parameter  sweeping. 

TRIMMED  FLIGHT 

In  finding  the  trim  point,  the  program  iterates  on  the  control  positions, 
fuselage  orientation,  rotor  attitude  and/or  engine  power  to  reach  desired  values 
of  the  rotor  flapping  moments,  forces  and  moments  on  the  aircraft  center  of 
gravity  and/or  engine  horsepower.  When  these  desired  values  are  achieved  the 
rotorcraft  is  trimmed. 

The  program  also  permits  the  calculation  of  two  steady  but  accelerated 
flight  conditions.  In  one,  the  rotorcraft  is  in  a  pushover  or  pullup  condition 
at  an  input  g-level.  In  the  other,  the  rotorcraft  is  in  a  banked  turn,  either 
level  or  spiral,  at  an  input  g-level.  In  these  cases,  the  desired  net  forces  and 
moments  on  the  rotorcraft  are  not  all  zero,  but  depend  on  the  user-requested 
g-level.  Either  an  unaccelerated  or  steady  accelerated  flight  condition  may  be 
used  as  the  starting  point  of  a  maneuver  simulation. 
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Figure  1.  Trim-Only  Operation 
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Figure  3.  Trim  or  Trim  and  Rotorcraft  Stability 
Analysis  Followed  by  Parameter  Sweep. 
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Figure  4.  Trim  Followed  by  Maneuver 
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Figure  5.  Trim  Followed  by  Maneuver  with  Rotorcraft  Stability  Analysis. 
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Figure  6.  First  Maneuver  in  Restart  Procedure 
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Figure  7.  Second  and  Lacer  Restart  Maneuvers. 
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A  complete  Input  data  deck  for  AGAP80  can  be  divided  Into  the  53  groups  or 
sets  of  cards  listed  In  table  I.  The  first  45  groups  form  the  basic  card  deck, 
which  is  used  for  trim-only  and  trim-and-rotorcraf t-stabillty-analysis-only 
program  operations.  The  remaining  eight  groups  are  only  included  In  the  deck 
when  a  maneuver  is  to  be  simulated. 

A  sample  input  data  deck  for  a  trim  followed  by  maneuver  is  presented  in 
Appendix  A.  Some  of  the  more  significant  Inputs  will  be  discussed  here,  but  for 
a  more  detailed  explanation  of  all  the  inputs  see  reference  1.  In  discussing  the 
input  data,  attention  will  be  directed  to  the  card  numbers  printed  in  columns 
73-80. 


NPART  Card 

In  order  to  compute  a  trimmed  flight  condition  the  NPART  card  (card  002) 
requires  a  "1"  in  column  1-2.  This  card  includes  the  primary  program  control 
variable,  NPART.  Permissible  values  of  NPART  on  this  card  are  1,2, 4, 5, 7,  and  8. 
The  value  of  NPART  specifies  the  type  of  operation  to  be  performed. 

1  “  Trim  only 

2  -  Trim  followed  by  maneuver  (maneuver  not  to  be  restarted) 

4  “  Trim  followed  by  maneuver  (maneuver  to  be  restarted) 

5  *  Maneuver  restart 

7  ■  Trim  followed  by  rotorcraft  stability  analysis 

8  -  Retrieve  maneuver  data  from  tape  for  analysis 

Program  Logic  Group 

The  Program  Logic  Group  (cards  006-013)  is  one  of  the  most  important  groups 
in  the  deck  because  it  contains  the  bulk  of  the  program  logic  controls.  It 
controls  which  groups  must  be  Included  in  the  deck  and  the  program  options  that 
will  be  used  in  the  computations.  This  group  consists  of  seven  cards  with  14 
integer  inputs  per  card.  Each  element  of  the  IPL  array  is  identified  as 
IPL(l-98).  Sections  2.3  and  4.3  of  reference  1  describe  in  brief  and  in  detail 
respectively  all  of  the  program  logic  inputs.  The  logic  inputs  control  the  three 
phases  of  the  C81  program: 

1)  Input  group  control  logic  -  controls  the  data  groups  which  must  be 
included  in  the  input  data. 

2)  Analysis  logic  -  controls  the  program  options  to  be  used  such  as 
unsteady  aerodynamics,  time-variant  rotor  analysis,  etc. 

3)  Output  control  logic  -  controls  the  data  to  be  output. 

The  logic  has  been  chosen  so  that  for  the  simplest  cases  most  Inputs  are 
zero.  In  general,  nonzero  inputs  activate  the  options  and/or  necessitate  inputs 
of  additional  data. 

Quasi-Static  vs.  Time-Variant  Trim  Analysis 

IPL  (49)  and  (50)  are  the  program  logic  Inputs  that  control  the  type  of  trim 
analysis  performed.  AGAP80  includes  three  types  of  trim  procedure: 

1)  Quasi-Static  (QS)  Trim 

2)  Quasi-Static  followed  by  Time-Variant  Rotor  (QS-TV)  Trim 

3)  Fully  Time-Variant  (FTV)  Trim 

12 
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TABLE  I  SEQUENTIAL  SUMMARY  OF  INPUT  GROUPS 
AGAP80  AGAP80 

SEQUENCE  SEQUENCE  REFERENCE 

NUMBER  OF  INPUT  CARD  CARD  SEQUENCE 

GROUP  TITLE _ ID  CARD  * _ NUMBERS  ** _ NUMBER  **** 


Deck  Identification  & 


Program  Flow  Control  Cards 

None 

00  -  04 

001  -  005 

Program  Logic  Group 

10 

10  -  17 

006  -  013 

Airfoil  Data  Table  Group 

Airfoil  Data  Table  Set  No. 

1 

21 

21/A, /B 1,/Cl, /D1 

014  -  295 

Airfoil  Data  Table  Set  No. 

2 

22 

22/A,/Bl,/Cl,/Dl 

296  -  506 

Airfoil  Data  Table  Set  No. 

3 

23 

23/A,/Bl,/Cl,/Dl 

N/A 

Airfoil  Data  Table  Set  No. 

4 

24 

24/A,/Bl,/Cl,/Dl 

N/A 

Airfoil  Data  Table  Set  No. 

5 

25 

25/A,/Bl,/Cl,/Dl 

N/A 

Airfoil  Data  Table  Set  No. 

6 

26 

26/A,/Bl,/Cl,/Dl 

N/A 

Airfoil  Data  Table  Set  No. 

7 

27 

27/A,/Bl,/Cl,/Dl 

N/A 

Airfoil  Data  Table  Set  No. 

8 

28 

28/A,/Bl,/Cl,/Dl 

N/A 

Airfoil  Data  Table  Set  No. 

9 

29 

29/A,/Bl,/Cl,/Dl 

N/A 

Airfoil  Data  Table  Set  No. 

10 

2A 

2A/A,/B1,/C1,/D1 

N/A 

Rotor  1  Group 

30 

30  -  3P 

507  -  524 

Rotor  1  Elastic  Pylon  Group 

40 

40  -  4V 

N/A 

Rotor  1  Elastic  Blade  Data  Group 

50 

50  -  5C 

525  -  765 

Rotor  2  Group 

60 

60  -  6P 

766  -  774 

Rotor  2  Elastic  Pylon  Group 

70 

70  -  7K 

N/A 

Rotor  2  Elastic  Blade  Data  Group 

80 

80  -  8C 

N/A 

Rotor  Aerodynamic  Group 

90 

90  -  9AE 

775  -  785 

Rotor  1  Rotor-Induced  Velocity 
Distribution  (RIVD)  Table 

100 

100  -  100/F 

786  -  1029 

Rotor  2  RIVD  Table 

110 

110  -  110/C 

N/A 

Ro  tor-Wake-at-Aerodynami c- 
Surface  (RWAS)  Table  Group 

RWAS  Table  No.  1 
RHAS  Table  No.  2 
RWAS  Table  No.  3 
RWAS  Table  No.  4 
RWAS  Table  No.  5 
RWAS  Table  No.  6 
RWAS  Table  No.  7 
RWAS  Table  No.  8 
RWAS  Table  No.  9 


*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 


*** 


*** 
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TABLE  I  (Continued) 


AGAP80 

AGAP80 

SEQUENCE 

SEQUENCE 

REFERENCE 

NUMBER  OF 

INPUT  CARD 

CARD  SEQUENCE 

GROUP  TITLE 

ID  CARD  * 

NUMBERS  ** *** 

number  **** 

RWAS  Table  No.  10 

*** 

*** 

N/A 

RWAS  Table  No.  11 

*** 

*** 

N/A 

RWAS  Table  No.  12 

*** 

*** 

N/A 

Basic  Fuselage  Group 

120 

120  - 

125 

1030  -  1032 

Fuselage  Aerodynamic  Group  or 
Fuselage  Aerodynamic  Table 

130 

130  - 

13C 

1033  -  1044 

Wing  Group 

140 

140  - 

14C 

1045  -  1055 

Stabilizing  Surface  Groups 
Stabilizing  Surface  No.  1 

150 

150  - 

15B 

1056  -  1065 

Stabilizing  Surface  No.  2 

160 

160  - 

16B 

1066  -  1077 

Stabilizing  Surface  No.  3 

170 

170  - 

17B 

1078  -  1089 

Stabilizing  Surface  No.  4 

180 

180  - 

18B 

N/A 

Jet  Group 

190 

190  - 

192 

N/A 

External  Store/Aerodynamic 

Brake  Group 

200 

200  - 

204 C 

N/A 

Rotor  Controls  Group 

210 

210  - 

21B 

1090  -  1094 

Iteration  Logic  Group 

220 

220  - 

22B 

1095  -  1098 

Flight  Constants  Group 

230 

230  - 

234 

1099  -  1103 

Bobwelght  Group 

240 

240  - 

241 

N/A 

Weapons  Group 

250 

250  - 

251 

N/A 

SCAS  Group 

260 

260  - 

264 

N/A 

Stability  Analysis  Times  Group 

270 

270  - 

272 

N/A 

Blade  Element  Data  Printout 

Times  Group 

280 

280  - 

282 

N/A 

Maneuver  Time  Card 

None 

291 

1104 

Maneuver  Specification  Cards 

None 

301  - 

30Z 

1105  -  1110 

Maneuver  Analysis  Cards 

None 

GDAP80 

Input 

1111  -  1117 

Cards 


*  "None"  indicates  that  the  group  does  not  have  an  identification  (ID)  card. 

**  C81  Program  A6AP80  version  sequence  card  numbers  for  input  data  deck  as 
described  in  Reference  1. 

***  No  specific  sequence  number  on  RWAS  Table  Cards. 

****  Refers  to  card  sequence  number  in  columns  73  -  80  of  sample  input  data 
deck  in  Appendix  A. 

"N/A"  indicates  that  the  group  is  not  included  in  the  sample  input  data  deck. 
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The  terms  Quasi-Static  and  Time-Variant  refer  to  the  primary  rotor  analyses  that 
are  used  in  the  procedure*  Both  analyses  can  include  the  effects  of  blade 
elasticity  (mode  shapes  other  than  the  rigid  body  mode  shape) ,  but  only  the 
Time-Variant  analysis  can  include  the  interaction  of  the  aerodynamic  loads  and 
blade  elasticity,  i.e.,  aeroelasticity •  In  the  QS  trim  procedure,  the 
Quasi-Static  rotor  analysis  is  used  for  both  rotor  systems.  In  the  QS-TV  trim 
procedure,  the  standard  QS  trim  is  performed  first.  The  rotor  modal  equations  of 
motion  are  then  numerically  integrated  for  a  prescribed  number  of  rotor 
revolutions  using  the  Time-Variant  rotor  analysis  with  the  control  positions  and 
fuselage  orientation  held  fixed  at  the  positions  determined  by  the  QS  trim.  It 
is  assumed  that  after  a  number  of  revolutions,  the  aeroelastic  effects  included 
in  the  Time-Variant  analysis  will  have  caused  the  rotor  either  to  stabilize  or 
diverge,  depending  on  the  basic  stability  of  the  rotor  system.  The  user  may 
select  that  such  Time-Variant  trims  be  computed  for  either  or  both  rotors.  If 
both  rotors  are  selected,  the  two  Time-Variant  trims  are  computed  independently 
of  each  other.  In  the  FTV  trim  procedure,  the  user  specifies  the  rotor  or  rotors 
which  use  the  Time-Variant  analysis.  In  doing  so,  any  reference  to  the  QS  trim 
procedure  is  deleted  for  the  specified  rotor (s).  A  rotor  which  does  not  use  the 
Time-Variant  analysis  uses  the  Quasi-Static  analysis.  Table  II  shows  the  type  of 
rotor  analysis  used  for  each  rotor  as  a  function  of  the  values  of  IPL(49)  and 
(50). 


Trim  Procedure 

The  basic  program  flow,  which  is  the  same  for  all  three  types  of  trim,  is 
shown  in  figure  9.  The  three  basic  blocks  in  the  figure  are  the  computation  of 
the  partial  derivative  matrix  (shown  in  figure  10) ,  the  Time-Variant  portion  of 
the  QS-TV  trim  (shown  in  figure  11),  and  the  rotor  force  and  moment  computations 
(shown  in  figure  12). 

Trimmed  Flight  Printout 

Two  types  of  printouts  are  possible  for  the  data  computed  in  the  last  trim 
iteration,  the  standard  trim  page  (figure  13)  and  the  optional  trim  page  (figure 
14).  The  standard  trim  page  follows  the  final  trim  iteration.  The  final  trim 
Iteration  occurs  either  when  all  forces  and  moment  imbalances  are  within  their 
respective  allowable  errors  or  after  the  iteration  limit  for  trim  has  been 
reached.  These  allowable  errors  and  the  iteration  limit  are  Inputs  of  the 
Iteration  Logic  Group. 

Printout  of  the  optional  trim  page  is  controlled  by  IPL(73)  (Program  Logic 
Group).  The  optional  trim  page  is  most  useful  for  presenting  data  from  a  wind 
tunnel  simulation.  The  optional  trim  page  contains  the  following  information: 

1)  problem  identification;  2)  a  parameter  listing  of  rotor  controls,  rotor 
parameters,  wind  tunnel  parameters,  and  program  options;  3)  rotor  forces  and 
moments  listed  in  both  wind  reference  and  shaft  reference  systems,  and  4)  a 
summary  of  the  beam,  chord,  and  torsional  rotor  loads  is  printed  if  rotor  blade 
elastic  mode  shapes  have  been  Included  in  the  analysis. 
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TABLE  II  ROTOR  ANALYSIS  USED  DURING  TRIM  AND  MANEUVER 


Inputs 


Analysis  Used 


IPLC49) 

IPL(50) 

Rotor 

In  Trim 

In  Maneuver 

0 

(Ignored) 

(1  (Main) 

QS 

QS 

12  (Tail) 

QS 

QS 

0 

J  1 

QS-TV 

TV 

f 

l  2 

QS 

QS 

1 

^  1 

/  1 

QS 

TV 

1 

X  2 

QS 

QS 

2 

/  1 

TV 

TV 

l  2 

QS 

QS 

0 

i  1 

QS 

QS 

[ 

t  2 

QS-TV 

TV 

2  < 

1 

i1 

QS 

QS 

l  2 

QS 

QS 

2 

i  1 

QS 

QS 

\  2 

TV 

TV 

0 

J  1 

QS-TV 

TV 

t  2 

QS-TV 

TV 

3 

1 

/  1 

QS 

TV 

\  2 

QS 

TV 

2 

(  1 

TV 

TV 

(  2 

TV 

TV 

QS  “  Quasi-Static  rotor  analysis 
TV  “  Time-Variant  rotor  analysis 

QS-TV  ■  Quasi-Static  trim  followed  by  a  Time-Variant  trim.  During  the 
time-variant  portion  of  this  type  trim,  only  the  flapping,  pylon,  and  mast  windup 
angles  of  the  time-variant  rotor  are  allowed  to  vary;  the  fuselage  and  control 
positions  are  held  fixed  at  the  values  determined  by  the  quasi-static  trim.  If 
both  rotors  are  time-variant,  they  are  trimmed  independently  of  each  other. 
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PDM  =  Partial  Derivative  Matrix 
QS-TV  =  Quasi-Static  Trim  Followed 
by  Time-Variant  Rotor  Trim 


Figure  9.  Flow  Chart  of  Trim  Procedure 
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PDM  *  Partial  Derivative  Matrix 


Figure  10.  Flow  Chart  of  Partial  Derivative 
Matrix  Computation. 
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f  START  OF  \ 
y  TV  TRIM  J 


- ] 

STARTING  FR 
NUMERICAL!.' 
ROTOR  MODAL 
MOTION  FOR  F 
ROTOR  REV< 

l _ 

DM  QS  TRIM, 

{  INTEGRATE 
EQUATIONS  OF 
IVE  COMPLETE 
DLUTIONS* 

1 

COMPUTE  TU 
ROTOR  FORCES 
FROM  LAST 

'IE-AVERAGE 

AND  MOMENTS 
REVOLUTION 

— 

WRITE  TRIM  P 
WITH  OUTPUT 
REFLECT  THE  R 
TIME-VARIAN 

l _ 

AGE  OR  PAGES 
UPDATED  TO 
ESULTS  OF  THE 

T  ROTOR  TRIM 

WRITE  BLADE  B 
SUMMARY  FROM 
DURING  LAST  RO 

ENDING  MOMENT 
DATA  COMPUTED 
TOR  REVOLUTION 

END  OF  "\ 
_  TV  TRIM  ) 


♦Control  positions  and  fuselage  degrees  of  freedom  held 
fixed  at  the  final  values  computed  during  the  QS  trim. 


Figure  11.  Flow  Chart  of  Time-Variant  Rotor  Trim. 
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«  Rotor-Induced  Velocity 
PDM  *  Partial  Derivative  Matrix 
KK  ■  Thrust/v^  Looping  Variable 


QS  *  Quasi-Static  Rotor  Analysis 
TV  *  Time-Variant  Rotor  Analysis 
NREV  *  Number  of  Rotor  Revolutions 
in  TV  Analysis 


Figure  12.  Flow  Chart  of  Rotor  Force  and  Moment 
Computations  During  Trim  Procedures. 
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Figure  14.  Optional  Trim  Page. 
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ROTORCRAFT  STABILITY  ANALYSIS 

The  rotorcraft  stability  analysis  operation  uses  a  trim  point  or  a  time 
point  during  a  maneuver  as  its  Initial  condition.  The  stability  derivatives  are 
computed  by  making  small  independent  perturbations  to  each  of  the  rotorcraft 
rigid-body  degrees  of  freedom  and  primary  flight  controls,  computing  the  forces 
and  moments  in  the  perturbed  flight  condition,  subtracting  the  values  of  the 
forces  and  moments  in  the  initial  condition  from  those  at  the  perturbed 
condition,  and  finally  dividing  the  differences  by  the  appropriate  perturbations. 
Using  these  derivatives,  the  roots  of  the  rotorcraft  equations  of  motion,  mode 
shapes  associated  with  the  roots,  and  transfer  functions  for  the  rotorcraft 
system  are  computed. 

Input 

The  input  cards  required  are  the  same  as  for  trim  only  (NPART-1).  In  order 
to  compute  a  trimmed  flight  condition  followed  by  a  stability  analysis,  an 
NPART«7  card  is  required.  Note  that  if  the  Time-Variant  rotor  analysis  is 
activated  for  either  rotor,  a  rotorcraft  stability  analysis  cannot  be  performed. 

A  stability  analysis  should  not  be  performed  for  V  -  0.  That  flight  condition 
should  be  simulated  with  some  small,  nonzero,  airspeed  (typically,  0.001  knot). 

Output  of  Rotorcraft  Stability  Analysis 

The  operation  of  the  rotorcraft  stability  analysis  routine  (STAB)  depends  on 
the  numerical  evaluation  of  a  number  of  partial  derivatives  that  appear  in  the 
equations  of  motion  for  the  rotorcraft. 

Control  Partial  Derivative  Matrices 

Force  and  Movement  Derivatives 

Figure  15  presents  the  control  partial  derivative  matrix  as  printed  out  in 
three  ways  by  the  AGAP80  version  of  C81.  The  first  version  of  the  control 
partial  derivative  matrix  is  printed  with  units  of  pounds  per  inch  or  foot-pounds 
per  inch.  The  response  to  each  of  the  14  degrees  of  freedom  available  in  STAB  is 
evaluated  and  ratioed  to  be  the  response  to  a  1-inch  step  input  from  each  of  the 
four  controls.  The  rotor  flapping  angles  are  changed  to  reduce  the  rotor 
flapping  moments  to  less  than  the  allowable  error  if  the  rotor  degrees  of  freedom 
are  not  turned  on. 

Control  Derivatives  In  Terms  Of  Accelerations 

The  second  version  of  the  control  partial  derivative  matrix  contains  the 
same  information  as  the  first.  In  this  matrix,  the  force  and  moment  derivatives 
have  been  divided  by  the  appropriate  masses  or  inertias  to  give  the  units  of 
linear  or  angular  acceleration  per  inch  of  control.  These  numbers  may  be  thought 
of  as  the  accelerations  at  the  Instant  immediately  after  a  step  input  from  the 
controls. 

Conventional  Fixed-Wing  Nondimensional  Derivatives 

If  the  rotorcraft  does  not  have  a  wing  or  if  the  airspeed  is  less  than  1.0 
knots,  this  matrix  is  not  printed.  Refer  to  reference  2  for  the 
nondlmenslonalizing  factors  and  for  interpretation  of  the  first  six  rows  of  the 
third  matrix.  No  attempt  will  be  made  to  Interpret  or  explain  the  last  four  rows 
of  this  matrix  because  conventional  fixed-wing  concepts  do  not  apply  to 
helicopter  rotors  and  pylons. 

Partial  Derivatives  For  Rotorcraft  Stability  Analysis  Degrees  Of  Freedom 

The  next  pages  of  output  contain  detailed  information  used  for  the 
calculation  of  the  partial  derivatives  for  each  degree  of  freedom  that  is 
activated  in  STAB.  The  partial  derivatives  are  evaluated  in  the  same  order  in 
which  the  variables  are  listed  in  the  following  paragraph.  See  figure  16. 
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Figure  15.  Control  Partial  Derivative  Matrix  from  STAB 
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Figure  16,  Example  of  Partial  Derivatives  for 
STAB  Degrees  of  Freedom. 
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Rotorcraft  Stability  Analysis  Degrees  Of  Freedom 

At  the  top  of  each  partial  derivative  page  (figure  16)  is  a  list  of  the 
current  value  of  each  of  the  possible  degrees  of  freedom.  All  FUS  (fuselage) 
parameters  are  in  the  body  reference  system  and  all  M.R.  and  T.R.  (rotor) 
parameters  are  in  the  appropriate  shaft  reference  system.  By  a  comparison  of  two 
successive  pages,  it  is  possible  to  tell  which  variable  is  being  perturbed  and  by 
how  much. 

The  30  variables  which  may  be  perturbed  are  perturbed  in  the  following 
order: 

FUS.  U  ■  Velocity  in  the  X  direction  (FT/SEC) 

FUS.  W  ■  Velocity  in  the  Z  direction  (FT/SEC) 

FUS.  Q  =  Pitch  Rate  (DEG/SEC) 

FUS.  V  *  Velocity  in  the  Y  direction  (FT/SEC) 

FUS.  P  =  Roll  rate  (DEC/SEC) 

FUS.  R  =  Yaw  rate  (DEG/SEC) 

M.R.  PYLON  MODE  1  RATE  -  (DEG/SEC) 

M.R.  PYLON  MODE  2  RATE  *  (DEG/SEC) 

M.R.  PYLON  MODE  3  RATE  -  (DEG/ SEC) 

M.R.  PYLON  MODE  4  RATE  -  (DEG/ SEC) 

T.R.  PYLON  MODE  1  RATE  -  (DEG/ SEC) 

T.R.  PYLON  MODE  2  RATE  -  (DEG/ SEC) 

T.R.  PYLON  MODE  3  RATE  -  (DEG/ SEC) 

T.R.  PYLON  MODE  4  RATE  -  (DEG/ SEC) 

M.R.  F/A  FLAP  RATE  -  (DEG/ SEC) 

M.R.  LAT  FLAP  RATE  -  (DEG/ SEC) 

T.R.  F/A  FLAP  RATE  -  (DEG/SEC) 

T.R.  LAT  FLAP  RATE  -  (DEG/ SEC) 

M.R.  PYLON  MODE  1  DISP  -  (DEG) 

M.R.  PYLON  MODE  2  DISP  -  (DEG) 

M.R.  PYLON  MODE  3  DISP  -  (DEG) 

M.R.  PYLON  MODE  4  DISP  ■  (DEG) 
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T.R.  PYLON  MODE  1  DISP  -  (DEG) 

T.R.  PYLON  MODE  2  DISP  -  (DEG) 

T.R.  PYLON  MODE  3  DISP  -  (DEG) 

T.R.  PYLON  MODE  4  DISP  -  (DEG) 

M.R.  F/A  FLAP  DISP  -  (DEG) 

M.R.  LAT  FLAP  DISP  -  (DEG) 

T.R.  F/A  FLAP  DISP  -  (DEG) 

T.R.  LAT  FLAP  DISP  -  (DEG) 

Each  partial  derivative  page  (figure  16)  gives  rotor  performance  data,  a 
force  and  moment  summary,  and  a  delta  force  and  moment  summary.  The  forces  and 
moments  printed  here  are  computed  after  the  small  Increment  In  the  pertinent 
variable  has  been  made.  All  data  are  In  the  body  reference  system.  The  delta 
force  and  moment  summary  presents  the  changes  In  the  force  and  moment 
contributions  In  exactly  the  same  format  as  the  full  force  and  moment  summary. 
Each  number  In  this  data  block  Is  obtained  by  taking  the  corresponding  value  from 
the  force  and  moment  summary  Immediately  above,  less  the  corresponding  value  at 
the  trim  condition  or  at  the  current  maneuver  time  point. 

Rotorcraft  Stability  Partial  Derivative  Matrices 

Total  Partial  Derivative  Matrix 

A  summary  of  the  partial  derivatives  computed  from  the  data  on  the  previous 
pages  Is  then  printed  (figure  17).  Each  row  gives  the  partial  derivatives  of 
some  total  aircraft  force  or  moment,  as  labeled,  with  respect  to  the  perturbation 
variables  used. 

Rotor  Partial  Derivative  Matrix 

A  summary  of  the  rotor  partial  derivatives  computed  from  the  data  on  the 
previous  pages  is  printed  on  this  page  as  well  (figure  17).  Each  row  gives  the 
partial  derivatives  of  some  force,  moment,  or  flapping  angle,  as  labeled,  with 
respect  to  the  linear  and  angular  velocities  U,  V,  W,  P,  Q,  and  R. 

Mass,  Damping,  and  Stiffness  Matrices 

The  mass,  damping,  and  stiffness  matrices  (figure  18)  which  are  used  to 
calculate  the  rotorcraft  stability  characteristics  are  printed  next.  Refer  to 
Volume  I  of  reference  3  for  the  analytical  background  of  these  three  matrices. 

Stick-Fixed  Rotorcraft  Stability  Results 

The  rotorcraft  characteristic  equation,  with  controls  fixed,  is  solved  for 
its  complex  roots  and  associated  response  modes.  The  results  are  presented  in 
several  ways  (figure  18). 
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Figure  17.  Rotor  and  Total  Partial  Derivative  Matrices 
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Figure  18.  Stability  Matrices  and  Stick-Fixed  Stability  Results 
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Figure  18.  Concluded 
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The  real  and  Imaginary  parts  of  the  roots  of  the  rotorcraft  characteristic 
equation  are  printed  under  the  headings  REAL  and  IMAG.  The  units  are  radians  per 
second.  The  roots  may  be  used  to  form  the  denominator,  d(s),  of  the  frequency 
response  polynomial. 


d(s)  -  n  (s  -  REAL  +  IMAG  *j>  (  s  -  REAL  -  IMAG  *j) 
i-1  i  i  i  i 


where 


s  «  Laplace  operator 

II  *  Continued  product  notation 

j  *  v/rT~ 

n  *  Number  of  roots  printed 
i  -  Sequence  number  of  root  in  printout 


Note  that  in  the  case  of  complex  conjugate  parts  of  roots,  only  the  root 
with  the  positive  imaginary  part  is  printed. 


Each  root  or  pair  of  roots  generates  the  terms  in  one  factor  of  the 
denominator  of  the  Laplace  transfer  function,  D(s). 


where 


D(s) 

» 

TAU 

i 

DAMP 


n  2 

II  (TAU  *  s  +  DAMP  *  s  +  1) 
i-1  i  i 


2  2 

1/(REAL  +  IMAG  )  =*  1/  w 

i  i  i 


2  2 

-2  *  REAL  /  (REAL  +  IMAG  )  «  -2  f  / w 

i  i  i  i  n  t 


and  II  ,  n,  and  i  are  as  defined  above. 


For  the  oscillatory  roots  of  the  rotorcraft  characteristic  equation,  the 
PERIOD  of  the  damped  oscillation  is  given  in  seconds.  For  the  roots  with  a  zero 
imaginary  part,  the  period  Is  a  meaningless  concept,  so  a  zero  appears  in  the 
output. 

The  column  headed  TIME  TO  HALF-DBL  depends  on  the  value  of  the  real  part  of 
the  root.  If  the  real  part  is  negative,  the  time  to  half  amplitude,  in  seconds, 
is  printed.  If  the  real  part  Is  positive,  the  time  to  double  amplitude,  in 
seconds,  is  printed.  The  column  headed  CYCLES  TO  HALF-DBL  contains  the  number  of 
cycles  to  half  or  double  amplitude  based  on  the  damped  natural  frequency  (  to  » 
IMAG)  for  the  oscillatory  roots.  A  zero  is  printed  for  aperiodic  roots. 

The  UNDAMPED  NATURAL  FREQ,  wn,  is  based  on  the  absolute  value  of  the 
complex  root. 

u«  n  -  y  REAL2  + IMAG2 

Thus,  wn  is  defined  even  for  an  aperiodic  root.  The  calculated  value  of 
u>  is  given  both  in  radians  per  second  and  cycles  per  second.  The  DAMPING 
RATIO,  f,  in  combination  with  the  undamped  natural  frequency,  completely 
describes  the  root. 


f  -  REAL  /  oj  ;  REAL  »  f  to 

n  n 

For  a  stable  aperiodic  root,  the  damping  ratio  is  I.  For  an  unstable  aperiodic 
root,  the  damping  ratio  is  -I. 
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In  Che  scabiliCy  mode  shape  printout  (figure  19),  each  column  represents  one 
mode*  The  first  column  on  the  left  is  associated  with  the  first  root  printed, 
the  second  with  the  second  root,  and  so  forth.  Each  component  of  a  mode  shape 
has  a  relative  magnitude  (MAGN)  and  a  phase  angle  (PHASE),  The  normal  printout 
provides  for  eight  columns  (mode  shapes  of  roots).  If  more  than  eight  roots  are 
computed,  the  additional  roots  are  printed  in  the  same  format  below  the  first 
set.  Columns  after  the  last  root  are  set  to  zero. 

The  mode  shapes  associated  with  the  rotorcraft  characteristic  roots  are 
first  printed  as  normalized  with  respect  to  THETA,  then  as  normalized  with 
respect  to  PHI,  and  lastly  as  normalized  with  respect  to  the  largest  variable  for 
a  given  root.  In  all  three  sets  of  normalized  mode  shapes,  the  normalizing 
variable  always  has  a  magnitude  of  1,000  (nondlmensional)  and  a  phase  angle  of 
0.0  degrees.  The  reader  should  see  reference  1,  section  6  for  a  description  of 
the  variables  used  for  the  mode  shapes. 

Transfer  Function  Numerator 

Following  the  mode  shapes,  the  numerators  of  the  transfer  functions  (figure 
20)  for  aircraft  response  and/or  flapping  angles  as  specified  by  IPL(93)  are 
printed.  Table  III  gives  the  value  of  IPL(93)  required  to  print  the  numerators 
of  the  transfer  functions  computed  by  STAB.  For  each  of  the  numerators  printed, 
the  value  labeled  GAIN  is  the  constant  term  in  the  frequency  response  polynomial; 
STATIC  GAIN  is  the  gain  term  to  be  used  in  the  Laplace  transfer  function. 

The  complex  roots  of  the  frequency  response  polynominal  are  printed  in  pairs 
of  columns  labeled  REAL  and  IMAG.  Below  the  real  and  imaginary  roots  are  the 
corresponding  values  in  the  numerator  of  the  Laplace  transfer  function,  TAU  and 
DAMP.  The  numerator  of  the  Laplace  transfer  function  N(s),  may  be  written  as 
follows: 

N(s)  -  STATIC  GAIN  *  (TAU  *  s2  +  DAMP  *  s  +  1) 

k-1  k  k 

The  STATIC  GAIN  is  the  ratio  N(s)/D(s)  evaluated  for  s  -  0. 

The  frequency  response  polynomial  as 

n(s)  -  (GAIN)  *  n  I(s  -  REAL  +  IMAG  *  1) 
k-1  k  k 

where  k  -  sequence  number  of  root 

m  -  total  number  of  roots  printed 

Zero  roots  are  not  printed  for  either  the  stick  fixed  or  control  input 
solutions,  so  the  final  order  of  the  transfer  function  generated  as  described 
above  may  be  incorrect.  In  this  case  usually  one  more  "s"  in  the  denominator 
will  correct  the  situation.  The  need  for  this  correction  may  be  found  by 
inspecting  the  numerator  and  denominator  polynomials.  The  transfer  function  is 
correct  when  the  highest  power  of  'V  for  the  denominator  is  2  larger  than  that 
for  the  numerator. 

Frequency  Response 

The  frequency  response  of  the  transfer  functions  (figure  21)  is  tabulated 
following  the  transfer  function  numerator  printout.  The  data  listed  are  the 
frequency  in  hertz  and  radians  per  second,  the  gain  in  the  decibel  equivalent  of 
a  magnitude  in  degrees  per  inch  of  control,  and  the  phase  in  degrees.  The  range 
of  frequencies  is  0.01  to  100  radians  per  second.  Construction  of  a  Bode  plot 
for  each  transfer  function  is  greatly  simplified  with  these  data. 


(s  -  REAL  -  IMAG  *  j)] 
k  k 


only  if  IMAG  4  0 


Figure  19.  Mode  Shapes  of  Stability  Result 


NADC-812 90-60 


O  O  O 

a  a  o  o 

•  ••I 


Urn  U  *U  U  So  U  U  U  JO  o 
*5  3t  So'  So  It  2-  «!  U  ft  2 


►  •  90  M« 
in  in 

I  I  < 


•a  HO  ««  n«  NOt  ««  »9  I« 

«•  I«  »m  •-  *9  *7  «0  «  ' 

I  I  ‘ 


Y  Y  Y  Y 

ai  «i  w  mi 


o«  mi  «<•  ♦«  5*  22  a  *2  St 

■to  •  *  •  9  •  s  •  *  •  #4  OM  9  **2  ?  ' 

09  «•  MO  *9  Nfi  5 •  O  ♦* 

--  «(•  ««  ««  *<-  «•  *j  “3  *3 


YYiYYYYY 


So  No  2*  So  8-  No  So  o  •• 

gn  g«  gg  gg  04  M  *  MM  M  *  O  f4 

ot  ot  o  •  O  •  Mom  Mo»  rn  •  MO  O  *  • 

Sn  ~M  So  So  0  .  -n  •««  -♦  *o  Of 

13  •«  •-  ••  •-  •-  •* 


Y  Y  Y  Y  Y  Y 


So 

a. 

S3 

••9 

5i 

ooNI 

09 

P 

P 

81 

S3 

-9 

OO 

oo 

09 

N4 

O  « 

P 

90 

O  * 

to 

•  • 

•T 

1i 

mO 

•  oo 

1 

OO 
•  oa 

18 

22 

•  o 

o  oo 

22 

1 

•Y 

m 

M 

m 

0* 

Y 

Y 

Y 

Y 

Y 

So 

So 

to 

MO 

9 

09 

So 

OMMO 

»o 

?! 

Si 

is 

99 

MO 

8 

Si 

*t 

O  • 

si 

Si 

IS 

•  • 

•  M 

►  O 

wo 

•  Oto 

22 

22 

MN 
•  9 

23 

IS 

12 

Y 

HI 

si  sc 


Y  Y  Y 
o«  ««  "o  So  No 

o  •  n  •  mm  n  •  n* 


51  n£  •  •  «o  o  9  •  2  •  5*m  2  •  2  • 

oA  A  •  Of  M  •  O  OM  MO  M  •  90  22 

go  o«  nn  «g  •  o  mo  o«o  oo  - o  nm 

•  m  « M  •«■  Ml  «•  •  •"  •“  •*  •**  •** 


oo  on  oo  So  o  o«  So  Co  So  o« 

21  52  SI  31  I  8*.  31  Si  31  8°. 

22  1J  28  IS  22  2$  2*  Is  18  2C 


m  m  m  m 

1  Y  Y  Y 

No  So  So  So  o  -o 

Si  8<5  8*.  81  8  Si 

«-  -o  o-  oo  >o  no 

•r  •■*  *f  *r  "•  *r 


«»«*•» 

Y  Y  Y  Y 

•l  «  u 
-•o  OO  «l»  oo  «► 

82  SC  31  Si  Si 

o  no  o-  o«  oo  oo 

•  *T  *"  *f  *“ 


;  :  o  c  :  ;  : 

{  i  j  ;  ;  1  J  1  ^ 

«  5  *  *  *  2  .*  «*  i 

111* 


THETA  /  LONG  CYC  PHI  /  LAY  CVC  PSI  /  PEDAL 


NADC-81290-60 


■Nfl 

??? 


N44 

P 

N0pc»O  4  •cOi* 

Ndod 

N 

0>0<VI  —  "KMJlO 

*3*8 

A 

iNWirtNNOdW 

PC 

a  •  •  •* 

» 

I  ifl  ifl  (M  83  O  f»  P  n 

IA 

i  n«i»o*««4o 

c 

«NP  i«N4«NO 

e 

•••••' 

•  *>*«n  • 

• 

Q  *  ID  ••••»• 

e 

1  1 

N 


N 

■  nBki 

H 

? 

z 

???? 

W 

< 

W  WWW 

Z 

JlA«)APN<Oe» 

0 

c«o 

pPQiA 

PC 

Swh" 

«»»♦ 

< 

weeaeeepN 

u 

<«9" 

NONN 

9 

SP«C«<0  •  •  • 

PC 

1-  l$«l 

e»o« 

PNC  •  iNdIDS 

p 

*  iso 

•  •  IDJIDDM  • 

< 

PC  PC  •  « 

•  •  •  • 

■«Nn 

e 

oooo 

i 

Till 

(»«  SCO 

Ul 

0 

N«M«* 

pSSSW 

OCP  PN- 

N 

pone 

M 

« 

*  a*2  *%• 

» 

InSSn 

888*3 

a 

» 

ions  teepee 

PC 

<n«  in 

«e»ns 

O 

p 

« 

* 

• 

1 

«  *  •  in  t  t<*en 

• 

• 

O-  •*  • 

< 

a 

NM 

a 

«mw« 

w 

a 

?? 

z 

??? 

Will 

< 

www 

z 

z 

jpoecNOttcp 

<NnenMNpne 

_  a«e 

PC 

3  UNO 

cee 

3 

< 

wc«*eeep««e 

u 

p  «9n 

N  *•  <M 

i0«fldrt  •  tin  • 

PC 

P  ixff 

«n« 

Z 

CCCC  tCP  »N 

p 

0«p«00~NK 

t  t  •  ta«>*p«t 

< 

•  PC  •  • 

1 

1  II  T  i  1  M 

p 

l/l 

e 

oeooa 

i 

i  i  i  i  T 

PC 

w 

iu  uuiiiJiy 

e 

N  « 

« 

eo 

N 

coppianonp 

i 

n  io 

« 

<*« 

p 

encp«e«np 

w 

u«  e 

PC 

♦  • 

<2 

in«*»in*ND 

n 

<P  o 

• 

•  e 

nt 

8 

inoo  «of“0 

PC  •  •  »(M  «PC  • 

38 

• 

<  teCYPCNDilM 

s 

PC 

a 

N 

me 

a 

? 

z 

PC 

? 

?? 

w 

< 

W 

w  w 

z 

joptieceopc 

a 

p  e  r 

•  f 

PC 

<nnpcepe«*e 

3  n  e  P 

on 

e 

wPPMneneec 

u 

<p  O  P 

PN 

o 

(f «ND  •  1  ••  • 

^p  e  e 

♦  c 

ion  #  #epp  »p 

K 

•OO-O-ONN 

? f7777??7 

$ 

lA 

•  •  • 

35 


Figure  20.  Numerator  of  Transfer  Functions 
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TABLE  III  VALUES  OF  IPL(93)  TO  PRINT  THE  NUMERATORS  OF  THE  TRANSFER  FUNCTIONS 


Denominator 

Of  Transfer 
Function 

Numerator  of  Transfer  Function 

n 

D 

9 

V 

<P 

M/R  Flapping 

Long. 

Lat . 

Collective 

2 

1 

1 

3 

3 

3 

4,-1 

4,-1 

Long.  Cyclic 

2 

1 

0,-1 

3 

3 

3 

4,-1 

4,-1 

3 

3 

3 

2 

0,-1 

1 

4,-1 

4,-1 

Pedal 

3 

3 

3 

2 

1 

0,-1 

4,-1 

4,-1 

<iuo/acc»  SAiMioat  pwkidic)  cainiobi  PMAtaioesi  uimoat  pmieiuai 
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Figure  21.  Frequency  Response  of  Transfer  Function 
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TRIM  SWEEPS 

The  parameter  sweep  operation  (figure  3)  may  be  used  to  simulate  the 
stacking  of  TRIM  and  TRIM-STAB  data  decks  for  a  given  rotorcraft.  Within  a  sweep 
deck,  the  user  specifies  by  input  data  those  cases  in  the  sweep  for  which  a 
rotorcraft  stability  analysis  is  and  is  not  to  be  performed.  The  parameters  most 
frequently  swept  include  airspeed,  gross  weight,  center-of-gravity  station-line, 
incidence  of  an  aerodynamic  surface,  atmospheric  conditions,  and  g-level. 
Generally,  only  one  parameter  is  changed  from  case  to  case  within  a  single  sweep 
deck.  However,  any  number  and  combination  of  inputs  except  some  program  logic 
switches  and  the  values  in  some  data  tables  may  be  swept.  The  assumption  is  made 
that  each  desired  trim  condition  bears  some  relationship  to  the  previous  one,  and 
that  the  previous  trim  point  is  a  good  starting  condition  for  finding  the  next 
trim  point.  For  example,  in  a  speed  sweep,  a  change  of  20  or  30  knots  is  the 
most  that  should  normally  be  used  between  40  and  150  knots.  Outside  of  this 
range,  the  maximum  change  should  not  exceed  10  knots. 

Following  a  data  set  for  trim  only  or  trim-and-rotorcraft-stability  analysis 
(NPART  ■  1  or  7),  the  parameter  sweep  option  may  be  exercised  by  including  an 
NPART  *  10  card  after  the  Flight  Constants  Group  (card  234  -  labeled  with 
sequence  number  1103).  An  NPART  ■  10  card  permits  the  changing  of  user-selected 
inputs  and  retrimming  the  configuration.  On  the  same  card,  a  new  input  variable 
is  defined  (NVARA)  which  specifies  whether  a  trim  or  trim-and-rotorcraft- 
stability  analysis  is  to  be  performed  on  this  case  in  the  parameter  sweep.  If 
NVARA  ■  0,  the  program  will  attempt  only  to  iterate  to  a  new  trim  condition 
(equivalent  to  NPART  *1);  if  NVARA  f  0,  the  program  will  attempt  to  trim  and,  if 
successful,  will  also  perform  a  rotorcraft  stability  analysis  (equivalent  to 
NPART  -  7) . 

The  data  set  for  NPART  ■  10  consists  of  the  following  cards: 

First  Card:  Card  01  NPART  card  with  NPART  *  10 

Subsequent  Card(s):  $Change  Changes  to  input  data  using  namelist  input 

An  NPART  *  10  data  set  may  be  followed  only  by  another  NPART  *  10  data  set, 
and  postprocessing  program  (GDAP80)  inputs.  Figure  22  shows  an  example  of  input 
data  required  to  perform  an  airspeed  sweep  for  trim  only.  Input  parameter 
XFC(l),  forward  velocity,  a  parameter  in  the  Flight  Constants  Group,  is  changed 
each  time  throughout  the  sweep  and  the  case  is  retrimmed.  The  program  assumes 
the  last  trim  point  is  a  good  starting  point  for  the  next  trim  case.  The  reader 
will  note  that  these  NPART  ■  10  data  sets  follow  immediately  behind  card  234  of 
the  Flight  Constants  Group  (that  card  is  labeled  with  a  sequence  number  00110300 
in  columns  73-80  of  the  sample  data  deck  in  Appendix  A.) 

MANEUVER  SIMULATION 

All  maneuvers  require  a  trim  point  prior  to  computing  the  time  history  of  a 
maneuver.  The  trim  point  is  used  to  supply  the  initial  conditions  to  a  system  of 
differential  equations  that  describe  the  behavior  of  the  rotorcraft  in  a 
maneuver.  Various  external  inputs,  or  forcing  functions,  may  be  applied,  such  as 
control  movements,  gusts,  store  drops,  and  wing  incidence  change  independent  of 
control  motion.  At  times  specified  by  input  data,  the  maneuver  can  be  suspended 
while  a  rotorcraft  or  rotor  aeroelastic  stability  analysis  is  performed.  The 
maneuver  is  then  resumed  as  if  no  interruption  had  occurred  and  continued  until 
it  reaches  either  the  next  time  point  to  do  a  stability  analysis  or  the  end  of 
the  maneuver. 
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A  maneuver  restart  operation  is  begun  just  like  an  ordinary  maneuver  using  a 
trim  condition  as  a  starting  point.  The  only  difference  is  that  the  time-history 
variables  and  many  intermediate  variables  are  saved  on  the  restart  magnetic  tape. 
Subsequent  maneuver  restarts  use  the  condition  at  one  of  the  saved  time  points  as 
the  initial  conditions,  and  so  do  not  require  a  trim  condition  or  the  complete 
data  set  defining  the  rotorcraft. 

Figures  4  through  7  illustrate  the  type  of  operations  that  can  be  performed 


that  are  related  to  maneuver  simulation. 

Maneuver  Simulation 

_ Operation  NPART 

a)  Trim  followed  by  maneuver  2 

b)  Trim  followed  by  maneuver  2 

with  rotorcraft  stability  analysis 

c)  Trim  followed  by  a  maneuver  in  which  4 

maneuver  time-point  data  are  stored 

for  a  subsequent  restart  of  the  maneuver 

d)  Maneuver  restart  5 


An  NPART  *  2  card  specifies  the  computation  of  a  trimmed  flight  condition 
followed  by  a  maneuver.  Subject  to  the  Program  Logic  Group  (IPL)  values,  a  data 
set  of  cards  02  through  291  (the  maneuver  time  card)  plus  at  least  one  301-type 
card  (maneuver  specification  card  or  commonly  called  J-card)  must  follow.  The 
maneuver  start  time  on  card  291  is  set  to  zero  regardless  of  the  input  value. 

The  postprocessing  program  (GDAP80)  inputs  follow  the  last  J-card. 

Whenever  a  maneuver  simulation  is  called  for  (NPART  *  2,  4  or  5)  an 
additional  input  parameter,  NPRINT,  which  specifies  the  frequency  of  printout  of 
maneuver  data,  is  to  be  supplied  in  columns  4-6  of  the  NPART  card.  The  program 
prints  data  showing  initial  conditions  for  the  maneuver  (maneuver  time  t  *  0)  and 
for  every  NPRINTth  time  point  thereafter.  A  blank  or  zero  input  is  reset  to 
unity. 

An  NPART  *  4  will  do  the  same  as  NPART  *  2,  except  that  the  maneuver  data 
will  be  stored  so  Chat  it  can  be  recalled  at  a  later  date  for  a  maneuver  restart 
(NPART  ■  5).  The  use  of  this  option  will  require  the  assistance  of  the  local 
programmer  to  set  up  the  restart  tape. 

For  a  maneuver  restart  case  (NPART  ■  5),  the  maneuver  time  card  has  the 
start  time  input  as  the  time  at  which  the  maneuver  is  to  be  restarted.  It  must 
be  greater  than  zero  and  less  than  the  last  point  of  the  maneuver  being 
restarted.  The  time  for  restart  need  not  be  identically  equal  to  a  previous  time 
point.  Then,  at  least  one  maneuver  specification  (J-card)  is  required.  The 
inputs  to  GDAP80  follow  these  inputs. 
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Maneuver  Time  Card  and  Maneuver  Specification  Cards 

The  following  information  is  provided  on  the  maneuver  time  card  (card  291) 
as  input  data: 

TCI(l)  Start  time  of  maneuver  (seconds) 

TCI(2)  First  time  or  azimuth  increment  (seconds  or  degrees) 

TC1(3)  Time  to  stop  using  first  increment  (seconds) 

TCI('0  Second  time  or  azimuth  increment  (seconds  or  degrees) 

TCI (5)  Time  to  stop  using  second  increment  and  return  to  first  increment 
(seconds) 

TCI (6)  Time  to  stop  the  maneuver  (seconds) 

TCI (2)  is  used  to  specify  the  first  base  value  of  the  time  increment  (At) 
between  the  calculation  of  maneuver  time  points*  The  At  computed  from  TCI(2) 
will  be  used  during  the  interval  of  TCI(l)  to  TCI(3)  seconds  of  maneuver  time. 

If  TCI(2)  <  1.0,  the  input  is  taken  to  be  At  in  seconds.  If  TCI(2)  1.0,  the 
input  is  taken  to  be  the  increment  in  rotor  1  azimuth  location  in  degrees  between 
time  points;  in  this  case,  the  time  increment  to  be  used  is  defined  as 

At  -  TCI(2)  /  6  O 

where  fil,  is  the  rotational  speed  of  rotor  1  in  units  of  rpm  and  At  is  in 
seconds. 

It  may  be  desirable  to  change  the  value  of  At  because  of  a  change  in  rotor 
speed.  For  this  case,  TCI(4)  can  be  used  to  specify  the  value  of  A t  to  be  used 

between  TCI(3)  and  TCI(5)  seconds  of  maneuver  time.  Like  TCI(2),  TCI(4)  may  be 

either  a  time  or  azimuth  increment.  It  is  not  necessary  that  TCI (2)  and  TCI (4) 
be  the  same  type  of  Increment;  e.g.,  one  may  be  time  and  the  other  azimuth. 

If  TCI(6),  the  time  to  stop  the  maneuver,  is  greater  than  TCI(5),  the 
program  then  uses  the  At  based  on  TCI(2)  between  TCI(5)  and  TCI(6)  seconds  of 
maneuver  time.  If  TCI (5)  Is  the  time  to  stop  the  maneuver,  as  well  as  the  time 
to  stop  using  the  At  based  on  TCI(4),  the  TCI(6)  input  may  be  zero  or  blank.  If 
a  second  time  Increment  is  not  desired,  then  TCI(4)  and  TCI(5)  should  be  input  as 
0.0.  In  this  case,  TCI (6)  will  be  ignored  and  TCI (3)  is  taken  as  the  time  to 
stop  the  maneuver. 

If  NPART  -  2,  4,  or  5,  one  maneuver  specification  card  (card  301  or  J-card) 
must  be  included  and  up  to  20  may  be  included.  All  have  the  same  format  (II,  14, 

5X,  6F10.0).  It  is  not  necessary  to  have  the  J  values  in  numerical  order,  and 

there  may  be  several  cards  with  the  same  value  of  J.  It  is  necessary  that  "next 
J"  (in  col  1)  be  blank  on  all  of  these  cards  except  the  last  one,  which  must  have 
some  symbol  in  the  first  column. 

Permissible  values  of  J  are  from  1  to  45.  The  type  of  variation  that  occurs 
for  each  value  of  J  is  given  in  the  following  list. 
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J  a  1  Movement  of  collective  stick 

J  ■  2  Movement  of  longitudinal  cyclic  stick 

J  *  3  Movement  of  lateral  cyclic  stick 

J  ■  4  Movement  of  pedal 

J  =  5  Inactive 

J  =  6  Folding  rotors  aft  after  tilting  forward  and  stopping 

J  *  7  Inactive 

J  =  8  Inactive 

J  *  9  A  vertical  ramp  gust;  ramp  length  may  be  zero 

J  *  10  A  vertical  sine-squared  gust 

J  *  11  A  horizontal  ramp  gust;  ramp  length  may  be  zero 

J  ■  12  A  horizontal  sine-squared  gust 

J  =*  13  A  change  in  engine  torque  supplied 

J  *  14  A  change  in  auxiliary  thrust  supplied 

J  *  15  Inactive 

J  *  16  Weapon  fire 

J  *  17  Change  of  longitudinal  mast  tilt  angle  and  of  RPM  on  both  rotors 
J  *  18  Rotor  brake 

J  *  19  Inactive 

J  *  20  Sinusoidal  movement  of  controls  or  mast 
J  =  21  Inactive 

J  “  22  Inactive 

J  ■  23  RPM-dependent  hub  springs 

J  *  24  SCAS  roll  channel 

J  ■  25  SCAS  pitch  channel 

J  ■  26  SCAS  yaw  channel 

J  *  27  Folding  rotors  horizontally  after  stop 
J  ■  28  RPM  dependent  flapping  stops 

J  ■  29  Connecting  and  disconnecting  helicopter  controls 
J  ■  30  Rotor  moment  balancing  mechanism 
J  -  31  Changing  NPRINT 

J  *  32  Simplified  automatic  pilot  simulatior 

J  ■  33  Inactive 

J  *  34  Deployment  of  an  aerodynamic  brake 
J  ■  35  Dropping  an  external  store 

J  ■  36  Changing  incidence  or  control  surface  deflection  angles  of 
aerodynamic  surfaces 
J  ■  37  A  trailing  vortex  system 

J  *  38  Inactive 

J  *  39  Inactive 

J  ■  40  Inactive 

J  *  41  Roll  rate  input  to  autopilot  (P- Tracker) 

J  *  42  Pitch  rate  input  to  autopilot  (Q-Tracker) 

J  *  43  Yaw  rate  input  to  autopilot  (R-Tracker) 

J  ■  44  Normal  load  factor  input  to  autopilot  (G-Tracker) 

J  *  45  Rate-of-climb  input  to  autopilot  (RC-Tracker) 


The  input  format  for  each  of  the  currently  available  J-cards  is  given  in 
section  4.29  of  reference  1.  Start  and  stop  times  refer  to  the  time  from  the 
start  of  maneuver  unless  otherwise  noted. 
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Figure  23  presents  a  sample  input  for  a  maneuver  that  consists  of 
longitudinal,  lateral  and  pedal  doublets  input  sequentially  at  maneuver  time 
t  *  0.0,  1.0,  and  2.0  seconds  respectively.  The  data  set  required  follows 
immediately  behind  the  flight  constants  group  and  in  this  sample  case  consists  of 
cards  labeled  sequentally  1104  through  1110.  The  maneuver  time  card  (labeled 
1104)  indicates  that  the  overall  maneuver  begins  at  time  t  **  0.0,  ends  at  time 
t  ■  7.5  seconds,  and  the  time  increment  (  At)  between  the  calculation  of  maneuver 
time  points  is  .02  seconds. 

For  this  sample  case  the  following  input  format  defines  the  information 
contained  in  the  six  maneuver  specification  cards  (labeled  1105  -  1110)  needed  to 
simulate  this  maneuver. 

Col  1  Blank  except  for  last  card  in  group 

Col  2-5  J,  *  2  longitudinal,  -  3  lateral,  *  4  pedal  movement 


Col  11-20  Start  time  for  input  rate  1  (sec) 

Col  21-30  Input  rate  1  (in. /sec) 

Col  31-40  Stop  time  for  input  rate  1  (sec) 

Col  41-50  Start  time  for  input  rate  2  (sec) 

Col  51-60  Input  rate  2  (in. /sec) 

Col  61-70  Stop  time  for  input  rate  2  (sec) 


For  normal  control  rigging,  positive  control  rates  correspond  to  up  collective, 
forward  longitudinal  cyclic,  right  lateral  cyclic  and  up  rotor  2  collective. 

Maneuver  Simulation  Output 

The  program  prints  the  contents  of  the  maneuver  time  card  (card  291)  and  all 

maneuver  control  cards  (all  301-type  cards,  the  J-cards)  before,  starting  the  trim 

procedure  (figure  24).  A  program-supplied  title  for  the  action  caused  by  each 
J-card  is  included  to  the  left  of  the  numerical  inputs  of  the  J-cards.  This 
serves  as  a  record  of  the  type  of  maneuver  specified  as  well  as  a  quick  way  to 
check  the  input  data. 

Prior  to  performing  the  maneuver,  the  program  will  calculate  and  print  out 
the  trim  iteration  pages,  a  standard  trim  page  and  an  optional  trim  page. 

It  is  possible  to  print  out  data  computed  during  a  maneuver  at  specified 
time  points.  The  value  of  NPRINT  on  CARD  01  specifies  that  data  is  to  be  printed 
each  NPRINTth  time  point.  For  the  sample  case  that  was  run,  NPRINT  *  5  and  At 
(time  increment  between  maneuver  time  point  calculations)  »  .02  seconds.  Then, 
for  this  case,  the  maneuver  time  point  data  is  printed  out  every  .1  second. 

The  format  of  and  data  on  the  maneuver-time-point  page  are  identical  to 

those  of  the  standard  trim  page  with  the  following  exceptions.  The  problem 
identification  data,  trim  condition  specification,  and  atmospheric  parameters  are 
omitted;  some  data  in  the  aerodynamic  surfaces  printout  are  changed;  and  some 
data  are  added  at  the  top  of  the  page  to  the  body  and  ground  reference 
parameters  and  to  the  rotor  parameter  printouts.  Figure  25  illustrates  a 
typical  maneuver-time-point  printout  page.  The  maneuver-time-point  page  is 
followed  by  a  force  and  moment  summary  in  both  body  axis  and  wind  axis  for  that 
time  point  (figure  26).  The  format  of  the  summary  is  identical  to  the  summary 
printed  during  trim  iterations. 


c 

o 


NADC-8 12 90-60 


•a 

u 

U  CO 

<u  o 
> 

3  <D 
0)  £ 
G  *H 


cO 

a 

O)  *H 

>  4* 

3  -H 
0)  cj 

c  a) 
cx 


sr- 


£ 


w  u 


o 

©  o  o 

e  o 

O 

3  O 

o 

o 

o  o 

o  o 

O 

o 

o 

o 

o 

o 

o 

o  o  o 

o 

o 

o 

o  o 

o 

o 

o  o 

o  o 

o 

o 

o  o  o 

o  o 

O 

O  O 

o 

o 

o  o 

o  o 

o 

o 

o 

o 

o 

o 

o 

3  0  0 

o 

o 

o 

o  o 

3 

o 

o  o 

o  o 

3 

■* 

JlfiN 

®  o 

OHIVifl 

-* 

in  9 

r»  ® 

O'  o 

*- 4 

CM 

ro 

•* 

in 

<osa 

O'  o 

—  CM  CO 

-4-  in 

>0  N  CO  O' 

O 

•*» 

i*-  r-  ►-  »»-  r»  ® 

CD  CD 

9 

9 

9  9 

9  9 

9  O' 

O' 

O' 

O' 

O' 

O' 

O'  o*  © 

O' 

o 

o 

o  o 

O 

o 

o  e 

o  o 

o 

o  o  o 

o  o 

o 

3  O 

o 

o 

o  o 

o  o 

o 

o 

o 

o 

o 

o 

o 

0  3  0  0 

— 

—4  »-4 

»“4 

^  *"4 

^4  ««4 

— 

o 

9  0  0 

o  o 

o 

o  a 

o 

o 

o  a 

O  O 

o 

o 

9 

O 

o 

o 

o 

3  3  0 

o 

3 

o 

O  3 

o 

3 

3  O 

^4  ^4 

o  o 

O 

o 

o  o  o 

o  o 

o 

e  o 

o 

o 

o  o 

o  o 

o 

o 

O 

o 

o 

o 

o 

O  3  O 

o 

O 

o 

3  O 

o  o 

o  o 

3  3 

9 

♦ 

9 

z 

CM 

CM 

o  9 

■o 

• 

CO 

•-•CM 

o 

o 

CM 

o 

o  o 

*—  • 

o 

o  o 

o 

CM 

o  o 

o  o 

o 

o 

o 

• 

• 

9 

5 

6 

5 

6 
5 

• 

•  • 

• 

• 

o  • 

•  • 

• 

• 

co  o  • 

*•4 

f- 

• 

•  • 

•  • 

• 

3  O 

-J 

Ui 

o 

o  o 

o 

•  o 

o  o 

o 

_lrll_ 

o 

•  cm  in 

i 

CM 

9 

N  P4 

r—  cm 

h- 

cr 

z 

a. 

o 

o 

M 

CJ 

»- 

o  o 

< 

•  • 

O 

<n 

_i 

CM 

o  o 

• 

O 

o 

_l  o 

o 

o  o 

o 

o 

o  o 

o  o 

o 

LU 

o 

O  *“4  • 

• 

in  o 

O 

O  • 

o  o 

3 

1  • 

o  • 

• 

•  • 

• 

• 

•  • 

•  • 

• 

cr 

• 

•  •  in 

ro 

CM  O' 

•  o 

m  • 

• 

•  9 

—  3 

-4 

3  3 

3 

3 

3  3 

3  3 

3 

X 

3 

o  M 

1 

1  CM 

i 

3  -4- 

3  n 

3 

o 

CJ 

i* 

a. 

O' 

in 

£ 

3 

9 

_l 

o 

n 

o 

CM 

in 

o 

• 

o 

cr  r- 

9  CO 

o 

CM 

o  o 

CM  O 

cs  o> 

o 

o  o 

o  o  o  n 

O  CM  O 

o 

O  O  • 

• 

in  • 

4  0  4 

o  -4 

O 

CM  • 

9 

• 

•  • 

• 

•  o  o 

•  CM 

• 

• 

•  •  in 

^4 

• 

•  • 

•  • 

• 

•  o 

cr  • 

o 

o  o 

o  o 

•  • 

o  « 

o 

•-4 

o 

co  —  ru 

I 

r»  » 

o 

—  —  CM  CM 

ro 

Ul 


O' 

ro 


CD 

in 


<c 

in 


cv« 


a  c 


9 

9  9 

9 

z 

3  •  • 

CM 

«  9 

CM 

X 

• 

O  o  a  o 

i  r-  co 

0 

0 

0 

• 

1 

0 

0  9  0  0 

r*-  co 

o 

cn 

o  cl  in  © 

•  cr  o  •  o  © 

r*4 

9 

^4 

O  44 

o 

•  o 

•* 

•  cn  in  •  •  • 

•  •  •  • 

O  sf 

• 

•  3  •  © 

in  <9  O  3  *9 

• 

• 

• 

• 

• 

• 

o  • 

• 

o  m  o  o  o 

o  o  o  o 

•  • 

o 

c  o  o  in 

CM  —  in  —  9 

o 

»4 

2 

2 

3 

a. 

cr 

in 

< 

3 

cr 

►- 

t- 

o 

z 

z 

cr 

in 

«c 

O  CM 

(9 

r»  i- 

1—  © 

n 

CM  9 

in 

N  —  w  • 

cn  • 

o 

O 

*> 

M  • 

o  ro 

in  9  9  9  •  X  O 

Z  0  —  0 

o 

o 

• 

• 

o 

o 

-4  in  CM 

o 

0 

R 

B 

6 

0 

0 

•  CM  -4  in 

CM  O 

o 

_l  • 

•  •  3  —  in  o 

•  o  •  •  in  in 

• 

•  o 

o 

• 

• 

0—3 

• 

•  0  —  9  •  • 

—  •  o  — 

«  • 

• 

o  O'  o  9  cm  j  •  in  m  o  in  —  cm 

•  o 

o 

^4 

^4 

in 

in 

•  •  • 

e 

e  x  1  •  e  o 

|  <-4  •  • 

•  o 

o 

qt— co—  i  o  —  cm  in  i  — 

v-4 

^4 

1 

l 

i 

i 

t— 

z 

»- 

X 

z 

o 

X 

X 

o 

►4 

19 

—1 

c 

1- 

►4 

< 

— 1 

O' 

cs  in 

O' 

(9 

®o  cr 

U. 

o 

^4  • 

o  in 

o 

•  •  ©  Ul  © 

O 

CM 

3  0  0  0 

o  i  r»  ®  o  o 

•  CM  o  o 

o  © 

o 

1  9  in  —  •  1—  •  CM 

•  m  •  in 

c 

o 

O' 

o 

o 

O' 

•  CM  • 

• 

•  X  9  9  •  • 

—  •  •  CM 

•  • 

• 

X  • 

—  —  o  —  o  •  in  _j  o  • 

• 

• 

• 

• 

• 

• 

• 

•  •  • 

o  o  <  ro  •  e  o 

1  —  O  • 

1  o 

o 

<  r- 

i  i  —  —  o  cm  o  r>  i 

o  o 

o 

^4 

—  CM 

CM 

->  o 

"3 


r-  » 

m  ♦ 


CM 


N-  —  ©  O 


IS  •*  ©  o  o 

® 


•  •  e  o  — 


o  o 

•  •  in 

CM  CM  • 
M  H  D 


CM  CM  ro  co  •* 


•  CM  • 

—  •  in 

•#  ©  CM 


•  o  o 
9  •  • 
CM  in  CM 
—  !*»  I 


44 


Figure  23.  Maneuver  Simulation  Input  Data 
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Figure  24.  Maneuver  Specification  Printed  as  Output 
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Figure  25.  Maneuver-Time-Point  Printout  Page 
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Figure  27  illustrates  for  the  sample  case  a  time  history  of  the  sequential 
doublet  control  inputs  and  the  resulting  aircraft  attitudes  and  rates  in  all 
three  axes.  This  time  history  was  produced  by  hand  plotting  the  data  from  each 
of  the  maneuver-time-point  printout  pages  that  were  computed  for  this  case.  It 
will  be  shown  later  in  the  description  of  the  postprocessing  program,  GDAP80, 
that  this  data  can  be  plotted  automatically  by  the  computer. 

RETRIEVING  MANEUVER  TIME-HISTORY  DATA  STORED  ON  MAGNETIC  TAPE 

A  small  portion  of  AGAP80  is  invoked  to  read  a  maneuver  data  tape  that  had 
been  created  during  a  previous  maneuver  simulation.  The  data  read  from  the  tape 
is  transferred  to  a  disk  file  for  subsequent  postprocessing  by  the  postprocessing 
program,  GDAP80. 

An  NPART  -  8  card  causes  data  stored  on  tape  to  be  loaded  on  the  plot  disk 
(see  figure  8).  The  local  programmer  should  be  consulted  prior  to  its  use.  The 
value  of  NVARA  on  this  card  must  not  be  equal  to  zero.  The  three  comment  cards 
(cards  02,  03,  and  04)  and  the  GDAP80  inputs  constitute  the  remainder  of  the 
NPART  -  8  deck.  Comment  card  no.  1  (card  02)  has  a  numeric  title  for  the  data 
set  for  identification  purposes.  It  is  printed  in  the  output  heading  and  is 
titled  IPSN.  Note  that  the  IPSN  in  columns  4-10  of  card  02  must  be  identical  to 
the  IPSN  used  in  the  run  which  created  the  tape. 

GDAP80  OPERATIONS 

During  the  course  of  running  a  trim  or  maneuver,  the  values  of  a  large 
number  of  time-history  variables  at  each  time  point  are  saved  on  the  plot  disk. 

At  the  conclusion  of  the  trim  or  maneuver,  postprocessing  operations  specified  in 
the  GDAP80  input  data  are  performed  on  these  variables.  Eight  postprocessing 
operations  can  be  performed  on  data  created  by  AGAP80,  as  outlined  in  figure  28. 

A  brief  explanation  of  each  of  these  postprocessing  capabilities  is  given  in 
the  following  sections.  At  present  (time-history  plots  and  storing  time-history 
data  on  tape)  are  planned  to  be  utilized  by  this  group  but  many  of  the  others  are 
considered  valuable  analysis  tools  for  rotary  wing  performance  and  structural 
analysis  personnel. 

TIME-HISTORY  PLOTS  (FIGURE  29) 

The  user  may  request  time-history  plots  of  any  of  the  variables  saved  during 
trim  or  maneuver  simulations.  Section  9  of  reference  1  provides  a  table  of  over 
2300  variables  whose  values  are  saved  at  each  time  point  during  a  maneuver 
simulation.  Code  numbers  are  used  to  identify  the  variables  to  be  plotted.  The 
code  number  for  each  variable  saved  is  given  in  table  27  of  reference  1.  These 
plots  may  be  output  on  the  printer,  a  CALC0MP  drum  plotter,  or  both.  Plots 
created  as  part  of  a  maneuver  restart  run  will  contain  data  for  the  entire 
maneuver. 

AEROELASTIC  STABILITY  ANALYSES  (FIGURES  30  AND  31) 

Program  GDAP80  contains  two  analyses  capable  of  identifying  the 
frequency  and  damping  of  a  perturbed  multi-variable  system.  The  user  can 
select  either  the  Moving  Block  Fast  Fourier  Transform  Analysis  (NPART  -  6)  or 
the  Prony  Analysis  (NPART  ■  13)  to  examine  the  nonsteady-state  response  of 
either  rotor.  See  sections  12.3  and  12.4  of  reference  3,  Volume  I  - 
Engineer's  Manual  for  an  understanding  of  these  tvo  analyses.  The 
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Figure  27.  Sample  Time-History  of 
Maneuver  Simulation. 
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Figure  28.  GDAP80  Postprocessing  Operations  -  Performed 
on  Time-History  Data  Created  by  AGAP80 
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Figure  29.  Plotting  Operation 
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|  Figure  30.  Moving  Block  Fast-Fourier  Transform  Operation. 
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Figure  31.  Prony  Stability  Analysis  Operation 
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frequency  and  damping  characteristics  of  airframe  motion  may  also  be  investigated 
using  these  analyses,  but  that  is  not  their  primary  function. 

STORING  TIME-HISTORY  DATA  ON  TAPE  (FIGURE  32) 

If  it  is  desired  to  perform  additional  postprocessing  of  the  saved 
variables,  they  may  be  transferred  from  the  plot  disk  to  a  magnetic  time-history 

tape  (NPART  -  8).  The  data  on  the  tape  may  then  be  reloaded  to  the  plot  disk  for 

further  use  at  a  later  time.  Note  that  time-history  data  which  has  been  stored 

on  tape  can  be  retrieved  with  an  AGAP80  deck  with  NPART  *  8. 

HARMONIC  ANALYSIS  (FIGURE  33) 

A  complete  harmonic  analysis  may  also  be  made  for  any  of  the  saved  variables 
(NPART  -  9).  A  Fast-Fourier-Transf orm  technique  is  used  to  examine  a  broad  range 
of  frequencies.  This  option  is  especially  useful  for  studying  rotor  bending 
moments  and  related  variables. 

VECTOR  ANALYSIS  (FIGURE  34) 

Frequently,  maneuvers  are  run  where  one  of  the  controls  or  the  longitudinal 
mast  tilt  angle  is  varied  sinusoidally.  In  this  case,  the  vector  analysis 
operation  (NPART  *  11)  can  be  very  useful.  This  analysis  uses  the 
least-squared-errors  technique  to  fit  the  saved  data  to  a  curve  of  the  form 

F^(t)  *•  A-j_SIN  (  ut  +  d>j[)  + 

Then,  any  amplitude  ratios,  Aj/Aj,  and  phase  angle  differences,  <t>±  -  <p*,  may  be 
computed.  Lastly,  linear  combinations  of  the  variables  may  be  derived  In  the 
following  form: 

F  (t)  -  C  F  (t)  +  D  F(t)  +  E4 

i  i  j  i  k  i 

CONTOUR  PLOTS  (FIGURE  35) 

Rotor  aerodynamic  quantities  can  be  tabulated  versus  radial  station  and 
azimuth  and  plotted  in  plane-polar  format  using  this  option  (NPART  -  12).  The 
data  are  plotted  assuming  that  the  blade  stations  are  equally  spaced  along  the 
radius.  The  tabulations  and  plots  are  particularly  useful  for  displaying  the 
rotor  aerodynamic  environment. 

CREATING  A  DATA  TRANSFER  FILE  (FIGURE  36) 

This  option  permits  the  user  to  transmit  C81-generated  data  to  a  temporary 
file  accessed  by  the  DATAMAP  File  Creation  program  in  order  to  add  the  data  to 
the  DATAMAP  Master  File  (NPART  ■  15).  The  user  can  then  use  DATAMAP  to 
postprocess  the  C81  data  and  compare  it  with  test  data  also  resident  upon  the 
Master  File.  (The  DATAMAP  programs  are  described  in  detail  in  reference  4). 

Note  that  this  option  is  described  here  to  demonstrate  a  potential  program 
capability.  This  installation  does  not  presently  have  the  DATAMAP  programs 
installed. 

USER'S  GUIDE  TO  THE  INPUT  AND  OUTPUT  FOR  GDAP80 

Data  generated  by  AGAP80  can  be  postprocessed  using  program  GDAP80,  which  is 
automatically  Invoked  following  an  AGAP80  run.  All  Inputs  to  GDAP80  must  follow 
all  inputs  to  AGAP80.  Data  to  be  postprocessed  must  have  been  generated  by  the 
AGAP80  run. 
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Figure  32.  Operation  for  Storing  Maneuver  Time-History  Data  on  Tape 
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Figure  33.  Harmonic  Analysis  Operation. 
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Figure  34.  Vector  Analysis  and  Data  Reduction  Operation 
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Figure  35.  Contour  Plot  Operation. 
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Figure  36.  Creation  of  a  Data  Transfer  File. 
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The  data  to  be  postprocessed  by  GDAP80  have  been  stored  in  one  or  more 
Postprocessing  Data  Blocks  (PDB)  as  they  were  created  by  AGAP80.  Separate  PDBs 
are  created  from  one  of  three  sources: 

1)  a  Quasi-Static  trim  for  which  IPL(79)  ^  0 

2)  a  Time-Variant  trim  of  either  rotor 

3)  a  maneuver 

Note  that  IPL(79)  controls  the  storing  of  certain  Quasi-Static  trim  rotor 
variables  for  tabulation  and  contour  plots.  Data  are  stored  in  one 
Postprocessing  Data  Block  for  each  rotor  during  Quasi-Static  trim  whenever  IPL 
(79)  +  0. 

Since  the  data  in  any  one  PDB  is  generally  completely  independant  of  the 
data  in  any  other  PDB,  the  postprocessing  instructions  for  a  particular  PDB  are 
input  to  GDAP80  in  a  unique  set.  Each  such  set  can  contain  instructions  to 
perform  any  of  the  following  operations,  with  certain  restrictions: 

1)  Plot  time  histories  of  selected  data  (NPART  =  3) 

2)  Perform  a  stability  analysis  of  time-history  data  using  the  Moving  Block 
Fast  Fourier  Transform  procedure  (NPART  ■  6) 

3)  Store  maneuver  time-history  data  on  magnetic  tape  for  future 
postprocessing  (NPART  *  8) 

4)  Perform  an  harmonic  analysis  of  time-history  data  (NPART  =  9) 

5)  Perform  a  vector  analysis  of  time-history  data  (NPART  51  11) 

6)  Tabulate  or  contour  plot  rotor  aerodynamic  data  (NPART  =  12) 

7)  Perform  a  stability  analysis  of  time-history  data  using  Prony's  Method 
(NPART  =»  13) 

8)  Create  a  Data  Transfer  File  (NPART  »  15) 

The  NPART  *  14  card  is  used  to  terminate  the  list  of  postprocessing 
instructions  pertaining  to  a  particular  PDB.  Upon  reading  an  NPART  *  14  card, 
GDAP80  indexes  the  AGAP80-generated  file  of  data  to  the  beginning  of  the  next  PDB 
and  prepares  to  execute  the  next  set  of  postprocessing  instructions. 

As  mentioned  earlier,  only  two  of  these  eight  postprocessing  operations, 
time-history  plots  and  storing  time-history  data  on  tape,  are  planned  at  present 
to  be  utilized  by  this  group.  The  input  format  and  output  for  these  two 
operations  is  discussed  in  the  following  two  sections. 

PLOTTING  OF  TIME- HI STORY  DATA 

The  input  format  for  plotting  of  time-history  data  is  illustrated  below. 
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Card  11 


Col  1-2  NPART  (must  equal  14  to  move  to  the  next  Postprocessing  Data 
Block) 

Card  21 

Col  2  NPART  (must  equal  3  for  plotting) 

Col  4-6  NPRINT  print  control 


Cards  22A,  22B. . .  (One  for  each  set  of  plots  desired  -  10  maximum) 


Col 

1 

KEYS 

Col 

2  - 

5 

KV1 , 

Col 

7  - 

10 

KV2 , 

Col 

12  - 

15 

KV3 , 

Col 

20 

KEY 

Col 

31  - 

40 

SCI , 

Col 

41  - 

50 

SC2, 

Col 

51  - 

60 

SC3 , 

(Blank  except  for  the  last  22-type  card) 

Code  for  variable  1 

Code  for  variable  2 

Code  for  variable  3 

(1  =  Plot  on  printer  only) 

Minimum  scale  for  KV1 
Minimum  scale  for  XV2 
Minimum  scale  for  KV3 


Section  9  of  reference  1  provides  a  table  of  the  code  numbers  to  be  used  for  KV1, 
KV2 ,  and  KV3 . 


Whenever  time-history  data  are  available,  the  20-series  cards  may  be  used  to 
plot  the  data.  This  procedure  is  an  option.  If  it  is  not  to  be  used,  simply 
omit  the  20-series  cards.  The  data  may  be  plotted  on  the  computer  printer  or  put 
on  tape  for  plotting  by  the  CALCOMP  plotter.  The  local  programmer  should  be 
consulted  for  the  proper  setup  for  jobs  that  write  a  tape  for  CALCOMP  plotting. 

Time-history  data  are  stored  after  all  time-variant  trim  cases  and  for 
maneuvers  and  may  be  plotted  by  inserting  20-type  cards  in  the  data  deck.  If  a 
maneuver  simulation  is  being  performed,  a  21-type  card  and  up  to  10  22-type  cards 
are  placed  immediately  after  the  maneuver  specification  cards.  See  figure  37, 
cards  labeled  with  sequence  numbers  1111  through  1117,  for  illustration. 


Card  21 


Column  2  must  contain  the  integer  3  to  call  the  plotting  routine.  NPRINT 
specifies  that  the  first  and  every  NPRINTth  data  point  following  are  to  be 
plotted.  If  NPRINT  =  0,  it  is  reset  to  unity. 

Cards  22A,  22B,  etc. 

The  first  column  of  all  but  the  last  22-type  card  must  be  blank.  The  first 
column  of  the  last  22-type  card  must  contain  a  character  (a  slash  is 
recommended) . 

One  22-type  card  is  required  for  each  plot.  A  maximum  of  10  of  these  cards 
is  permitted  after  each  CARD  21.  Each  plot  may  contain  one  to  three  variables. 
The  first  three  inputs  (KV1 ,  KV2 ,  KV3)  on  a  22-type  card  are  the  code  numbers  for 
the  variable(s)  to  be  plotted.  The  code  numbers  must  be  integers.  If  only  one 
variable  is  to  be  plotted.  The  code  numbers  must  be  in  columns  2-5;  if  only 
two  are  to  be  plotted,  only  columns  2-5  and  7-10  are  to  be  used. 

KEY  (column  20)  controls  where  the  plotting  is  done. 
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Figure  37.  GDAP80  "Input  Data  for  Plotting 
of  Time-History  Data. 
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KEY  »  0  for  CALCOMP  only 

KEY  =  1  for  printer  only 

KEY  *  2  for  both 

The  program  internally  computes  its  own  scales  for  plotting  each  variable 
based  on  the  maximum  and  minimum  values  of  the  variables  during  the  time-history 
and  internally  specified  minimum  scales.  The  internal  minimum  scale  may  be 
overriden  for  each  variable  with  the  last  three  inputs  on  the  22-type  card.  The 
minimum  scale  inputs  (SCI,  SC2,  SC3)  are  in  units  of  the  appropriate  variable  per 
inch  for  printer  plots  and  units  per  centimeter  for  CALCOMP  plots. 

If  the  user  wishes  to  plot  more  variables  than  permitted  on  the  22-type 
cards,  then  another  21-type  card  should  be  inserted  in  the  instruction  set, 
followed  by  up  to  10  more  22-type  cards. 

A  sample  time-history  printer  plot  that  was  generated  after  maneuver 
simulation  is  provided  in  figure  38.  The  plot  symbols  used  are  the  numbers  1,  2, 
and  4.  The  variables  corresponding  to  each  symbol  and  its  units  are  printed  as 
part  of  the  plot  heading.  If  two  or  all  three  of  the  curves  intersect  at  a 
single  point,  the  symbol  printed  is  the  sim  of  the  individual  symbols.  For 
example,  the  symbol  7  (=1+2+4)  means  that  all  three  curves  pass  through  the  point 
where  the  7  is  printed. 

The  lower  and  upper  limits  on  the  plot  scale  are  given  for  each  variable 
plotted.  The  scale  in  units  per  inch  is  also  given. 

The  user  is  cautioned  that  the  automatic  plot  scaling  procedure  may  expand 
small  variations  completely  out  of  proportion  to  their  true  importance.  Be 
certain  to  check  the  scales  on  all  plots. 

Time  is  the  independent  variable,  and  is  printed  along  the  left  edge  of  the 
plot,  defining  the  time  axis.  Maneuver  plots  will  always  start  at  t  =  0.0.  If 
the  time  increment  is  changed  at  some  point  during  a  maneuver,  there  will  be  a 
change  in  the  time  scale  at  this  point.  The  resulting  compression  or  expansion 
of  the  time  scale  may  cause  apparent  discontinuities  that  are  not  actually  in  the 
data.  The  user  should  check  the  time  scale  carefully. 

Each  plot  card,  CARD  22,  is  independent  of  all  other  plot  cards.  Thus,  if 
desired,  one  variable  may  be  plotted  on  more  than  one  plot. 

The  dots  printed  down  the  page  are  spaced  at  1-inch  intervals  to  make  it 
easier  to  read  the  plot  values  by  eye.  They  also  provide  reference  lines  to  help 
see  slower  variations  on  long  time  histories. 

The  plot  routine  can  store  the  values  of  all  plot  variables  for  a  maximun  of 
2000  time  points.  Should  the  user  specify  NPRINT  *  1  on  CARD  21  for  a 
particularly  long  maneuver,  the  program  will  keep  internally  doubling  NPRINT 
until  the  total  number  of  points  to  be  plotted  is  less  than  or  equal  to  2000. 

For  CALCOMP  plots,  the  names  of  the  variables  plotted  appear  at  the  top  of 
each  CALCOMP  page  along  with  their  respective  plot  symbols.  The  vertical  scales 
and  the  plots  themselves  are  identified  by  the  plot  symbols. 
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Figure  38.  Sample  Time-History  Plot 
after  Maneuver. 
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It  is  recommended  that  the  user  avoid  plotting  periodic  variables  of 
approximately  the  same  magnitude,  with  near-zero  phase  shifts,  on  the  same  plot, 
as  it  will  be  difficult  to  differentiate  between  the  traces.  This  applies  to 
printer  plots  as  well.  A  maximum  of  2500  time  points  may  be  plotted  on  CALCOMP 
PLOTS. 

STORING  TIME-HISTORY  DATA  ON  TAPE 

The  input  format  for  this  postprocessing  option  is  given  below. 

Card  41 


Note  that  this  card  may  be  included  only  for  the  Postprocessing  Data  Block 
resulting  from  a  maneuver. 

Col  2  NPART  (must  equal  8  for  tape  file  operations) 

Col  11-  15  NVARA  (must  be  blank  or  all  zeros  to  store  data) 

Note  also  that  time-history  data  which  has  been  stored  on  tape  can  be  retrieved 
with  an  AGAP80  deck  with  NPART  *  8. 

Following  a  maneuver  (NPART  *  2,  4,  or  5),  it  may  be  desirable  to  store  the 
time-history  on  tape  so  that  the  data  can  be  recalled  later  for  additional 
analysis  or  plotting.  Inputs  of  8  and  0  in  columns  2  and  15  respectively  will 
store  the  data.  However,  the  local  programmer  should  be  consulted  for  the  proper 
setup  of  the  job  before  attempting  to  use  this  option.  This  instruction  should 
only  be  used  for  the  PDB  resulting  from  the  maneuver  portion  of  a  run. 
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PRESENT  OPERATIONAL  PROBLEMS 

Two  program  operation  problems  were  identified  while  running  test  cases  to 
check  out  each  of  the  program  operational  capabilities*  The  two  problems 
encountered  were  1)  the  inability  to  perform  any  of  the  postprocessing  operations 
of  GDAP80,  and  2)  the  inability  to  compute  the  correct  eigenvalues  in  the 
rotorcraft  stability  analysis.  Both  of  these  problems  are  associated  with  the 
copy  of  the  AGAP80  version  that  is  installed  on  the  CDC  CYBER  175/CDC  6600 
co*nputer  system  at  NAVAIRDEVCEN  and  are  not  considered  universal  problems  with 
the  C81  program. 

POSTPROCESSING  OPERATIONS 

After  running  a  test  case  to  perform  a  maneuver  simulation,  it  was  decided 
to  check  out  the  postprocessing  capabilities  available  in  GDAP80,  namely,  the 
plotting  of  time-history  data.  Upon  submitting  the  test  case  with  all  the  proper 
required  inputs  to  plot  the  time-history  of  the  selected  parameters  (as  Indicated 
in  figure  37) ,  the  program  calculated  and  printed  out  all  the  maneuver-time-point 
digital  data  but  did  not  print  the  printer  plots  as  directed.  After 
investigating  the  problem  it  was  determined  that  the  CDC  conversion  copy  of  the 
AGAP80  version  of  C81  as  delivered  to  NAVAIRDEVCEN  did  not  have  the 
postprocessing  program  (GDAP80)  subroutines.  Further  checking  with  the  Army 
Structures  Laboratory  revealed  that  the  postprocessing  program  portion  of  C81  had 
not  been  converted  as  yet  for  usage  on  CDC  CYBER  computer  systems.  The 
conversion  of  this  part  of  the  program  is  presently  underway.  As  soon  as  the 
conversion  is  complete  and  validated  by  the  Army  Structures  Laboratory,  the 
NAVAIRDEVCEN  will  receive  a  new  copy  of  AGAP80  program  source  on  magnetic  tape 
which  will  Include  the  postprocessing  program  GDAP80.  It  is  anticipated  that 
this  work  will  be  completed  by  March  1982. 

EIGENVALUE  CALCULATION 

Several  test  cases  were  run  in  an  attempt  to  check  out  the  rotorcraft 
stability  analysis  capability  of  the  program.  Throughout  all  attempts  the 
program  would  successfully  calculate  and  print  out  the  stability  derivatives  for 
both  the  total  aircraft  and  the  rotors.  The  program  would  then  develop  and  print 
out  the  mass  matrix,  the  damping  matrix,  and  the  stiffness  matrix.  Up  to  this 
point  all  output  data  correlated  with  the  sample  rotorcraft  stability  analysis 
case  in  section  6.8  of  reference  1.  But,  when  the  controls-f ixed  roots  and 
transfer  function  denominator  were  printed  out,  it  was  noted  that  that  output 
data  had  to  be  Incorrect.  For  example,  the  printout  would  one  time  show  the 
following  12  complex  roots  of  the  rotorcraft  characteristic  equation: 


Root  # 

Real 

Imag 

1 

•10000E+13 

0. 

2 

•10000E+13 

0. 

3 

•10000E+13 

65536 

4 

.10000E+13 

65536 

5 

.  10000E+13 

65536 

6 

•10000E+13 

65536 

7 

.  10000E+13 

65536 

8 

•10000E+13 

65536 

9 

.  10000E+13 

65536 

10 

.10000E+13 

65536 

11 

•  IOOOOE+13 

65536 

12 

.10000E+13 

65536 
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and  the  program  would  stop  printing  at  that  point.  That  is,  it  would  not  go  on 
to  print  out  the  remaining  parts  of  the  stability  analysis,  namely,  the 
controls-f ixed  mode  shapes,  transfer  function  data,  and  frequency  response  data. 

Another  run  of  the  same  case  would  produce  different  results.  An  example  of 
different  results  is  the  following  4  complex  roots  of  the  rotorcraft 
characteristic  equation: 


Root  # 

Real 

Imag. 

1 

. 10000E+13 

0. 

2 

.10000E+13 

0. 

3 

. 10000E+13 

.  15749E-03 

4 

.10000E+13 

518.96 

and  then  the  program  would  continue  printing  out  extraneous  controls-fixed  mode 
shape  data,  transfer  function  data,  and  frequency  response  data  to  correspond 
with  these  4  strange  roots. 

There  have  been  several  discussions  about  this  problem  with  personnel  at  the 
Army  Structures  Laboratory  and  the  Applied  Technology  Laboratory,  Ft.  Eustis,  VA. 
This  problem  is  still  under  investigation  and  is  believed  to  be  associated  with 
the  conversion  of  the  code  for  usage  on  CDC  CYBER  computer  systems.  The  two 
subroutines,  ALLMAT  and  ALLVEC,  which  the  stability  analysis  uses  to  compute 
eigenvalues  and  eigenvectors  have  had  some  cards  changed  as  a  result  of  the 
conversion.  It  is  expected  that,  with  the  second  conversion  of  AGAP80  that  is 
currently  underway  by  the  Army  Structures  Laboratory  to  include  the 
postprocessing  program  GOAP80,  this  problem  will  be  resolved. 
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COST  OF  PROGRAM  OPERATION 

The  cost  of  running  this  program  on  the  CD C  CYBER  175/CDC  6600  computer 
system  at  NAVAIRDEVCEN  has  been  of  interest.  With  the  limited  number  of  cases 
that  have  been  run  thus  far  throughout  this  checkout  and  validation  stage,  a  cost 
breakdown  for  trim  analysis,  stability  analysis,  and  maneuver  simulation  is 
presented  in  the  following  sections. 

TRIMMED  FLIGHT  ANALYSIS 

The  cost  of  running  a  trimmed  flight  analysis  is  directly  related  to  the 
number  of  Iterations  required  to  converge  to  a  trimmed  flight  condition.  Two 
factors  determine  this  number  of  iterations.  The  program  is  set  up  to  accept 
inputs  of  allowable  errors  in  the  force  and  moment  summations  (XIT(15)  through 
XIT(21)  of  the  Iteration  Logic  Group).  Normally,  the  smaller  the  allowable 
errors,  the  larger  the  number  of  iterations  required  for  trim. 

Also,  the  initial  settings  of  the  control  positions,  the  Euler  angles,  the 
rotor  flapping  angles  and  rotor  thrust  inputs  have  a  direct  bearing  on  how  many 
iterations  are  required  to  converge  to  a  trim.  These  are  all  Inputs  in  the 
Flight  Constants  Group  (XFC(S)  through  XFC(2Q)).  The  closer  these  values  are  to 
the  final  trim  value,  the  faster  the  analysis  will  converge  on  trim.  In  a 
parameter  sweep,  the  previous  trim  point  is  the  starting  condition  for  finding 
the  next  trim  point.  So,  for  example,  in  a  speed  sweep,  a  change  of  20  or  30 
knots  is  the  most  that  should  normally  be  used  between  40  and  150  knots.  Outside 
of  this  range,  the  maximum  change  should  not  exceed  10  knots. 

A  typical  trimmed  flight  analysis  cost  about  $25  to  $30  and  used  about  20  to 
25  cp  seconds  execution  time.  Whereas,  a  parameter  sweep  of  five  airspeeds  for 
trim  cost  about  $95  to  $100  and  used  about  75  to  80  cp  seconds  of  execution 
time. 


I 


I 


The  program  includes  the  option  for  two  basic  types  of  rotor  analyses  for 
trim:  Quasi-Static  (QS)  and  Time-Variant  (TV).  The  cost  data  discussed  above  is 
for  a  QS  rotor  analysis  trim.  A  Time-Variant  (TV)  trim  is  much  more  expensive 
and  is  generally  used  to  simulate  an  aeroelastic  rotor  when  blade  load  data  are 
desired.  The  Quasi-Static  (QS)  rotor  analysis  is  best  suited  to  rotor  systems 
where  the  dominant  blade  motions  are  at  1/rev  and  where,  for  most  performance  and 
flying  qualities  investigations,  other  flapping  frequencies  can  be  neglected. 

The  major  limitation  of  the  Quasi-Static  analysis  is  that  only  the  static  and 
1/rev  component  of  the  blades  response  are  calculated.  Consequently,  accurate 
computation  of  blade  loads  is  not  possible.  For  a  more  detailed  explanation  of 
the  two  types  of  rotor  analyses,  the  reader  is  referred  to  Volume  I,  section  2.2 
of  reference  3* 


Two  Identical  trimmed  flight  analyses  were  run  with  the  only  difference 
being  the  requirement  for  one  to  be  a  QS  trim  followed  by  a  TV  trim  of  the  main 
rotor.  The  data  below  gives  an  indication  of  the  difference  in  cost  to  perform 
these  two  analyses. 

Case  #  Type  of  Trim  Execution  Time  r*  CP  Seconds  Cost 


$86.10 

$27.72 


[»] 


QS  +  TV 
QS  only 


68.427 

21.555 
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The  Time-Variant  portion  of  a  Quasi-Static  trim  followed  by  a  Time-Variant 
trim  is  in  essence  a  time-history  of  XIT(6)  rotor  revolutions  with  the  fuselage 
and  control  positions  locked.  (In  this  case  XIT(6)  *  7  rotor  revolutions).  For 
each  rotor  which  is  Time-Variant,  the  additional  run  time  for  the  Time-Variant 
trim  after  the  Quasi-Static  trim  will  be  about  the  same  as  the  time  for  a 
maneuver  of  XIT(6)  rotor  revolutions. 

ROTORCRAFT  STABILITY  ANALYSIS 

A  rotorcraft  stability  analysis  can  be  performed  following  a  Quasi-Static 
trim.  Uith  the  eigenvalue  calculation  problem  that  was  discussed  in  an  earlier 
section,  only  preliminary  cost  data  for  this  operation  can  be  discussed.  At 
least  five  cases  were  run  in  an  attempt  to  checkout  this  analysis  operation.  In 
all  cases  the  program  stopped  short  of  computing  the  correct  eigenvalues.  Some 
cases  computed  incorrect  eigenvalues  and  used  those  in  further  computations, 
while  other  cases  caused  the  program  to  terminate  after  calculating  and  printing 
out  the  extraneous  roots. 

Even  with  these  problems  in  the  rotorcraft  stability  analysis,  all 
indications  lead  to  the  assumption  that,  if  the  analysis  were  completed  with 
correct  computations  of  the  eigenvalues,  the  cost  of  operations  would  be 
approximately  the  same.  The  following  data  gives  an  indication  of  anticipated  cost 
of  operation  of  the  trim  followed  by  rotorcraft  stability  analysis. 

Execution  Time 


Case# 

Type  Of  Analysis 

CP  Seconds 

Cost 

Comments 

6 

TRIM+STAB.ANAL. ( 1) 

42.067 

$52.77 

Program  terminated  before 
analysis  completed 

7 

TRIM+STAB.ANAL. ( 1) 

41.951 

$52.58 

Program  terminated  before 
analysis  completed 

8 

TRIM+STAB.ANAL  .  ( 1 ) 

43.363 

$54.34 

Analysis  completed  with 
incorrect  eigenvalues 

9 

TRIM+STAB.ANAL. (2) 

23.203 

$29.94 

Analysis  completed  with 
incorrect  eigenvalues 

12 

TRIM+STAB.ANAL. (2) 

26.585 

$33.93 

Program  terminated  before 
analysis  completed 

15 

TRIM+STAB.ANAL. (2)  59.487 

(3  Case  Parameter  Sweep) 

$74.04 

Program  terminated  before 
last  analysis  completed 

(1) 

Coupled  fuselage  (one  6x6  matrix) 
flapping  degrees  of  freedom 

in  stability  analysis  matrix;  no  rotor 

(2)  Uncoupled  fuselage  (two  3x3  matrices)  in  stability  analysis  matrix;  both 
rotor  flapping  degrees  of  freedom 


MANEUVER  SIMULATION 

The  amount  of  computer  execution  time  required  to  perform  a  maneuver 
simulation  is  dependent  on  two  factors:  1)  the  time  Increment  (At)  between  the 
calculation  of  maneuver  time  points;  and  2)  the  total  length  of  maneuver  time. 
These  two  factors  can  be  lumped  into  one  by  the  simple  fact  that: 
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number  of  maneuver  time  points 


total  length  of  maneuver  time 


(At)  time  increment  between  maneuver  time  points 


Therefore,  the  cost  of  running  a  maneuver  simulation  is  directly  related  to  the 
number  of  maneuver  time  points. 


Section  4.28  of  reference  1  discusses  some  of  the  contraints  placed  on  the 
size  of  At.  These  requirements  generally  apply  to  Time-Variant  maneuvers. 


The  sample  maneuver  simulation  case  discussed  earlier  was  a  trim  followed  by 
maneuver,  both  of  which  used  a  Quasi-Static  rotor  analysis.  The  maneuver  begins 
at  time  t  -  0.0,  ends  at  time  t  ■  7.5  seconds,  and  the  time  Increment  (At) 
between  the  calculation  of  maneuver  time  points  is  .02  seconds. 


A  At  equal  to  .02  seconds  was  found  to  be  both  adequate  and  necessary  for 
this  aircraft  with  this  rotor  system  with  a  rotational  speed  of  324.1  RPM,  using 
the  Quasi-Static  rotor  analysis  for  maneuver.  A  At  of  .02  seconds  corresponds 
with  a  maneuver  time  point  calculation  every  38.9  degrees  of  main  rotor  azimuth. 

It  was  found  that  when  At  was  input  as  .05  seconds  (every  97.2  degrees  of  rotor 

azimuth)  the  numerical  integration  technique  vent  unstable.  A  At  of  .02  seconds 
or  less  (approximately  30  degrees  or  rotor  azimuth  or  less)  is  recommended  for 
Quasi-Static  maneuvers.  To  insure  stability  of  the  numerical  technique  during  a 

Time-Variant  maneuver,  the  azimuth  increment  should  always  be  less  than  or  equal 

to  15  degrees. 


The  following  data  illustrates  the  cost  of  performing  maneuver  simulations 
with  this  program.  It  must  be  mentioned  that  the  cost  figures  are  for  a  digital 
listing  of  the  maneuver  time  point  data  and  do  not  represent  any  postprocessing 
of  maneuver  data,  such  as  time  history  plots,  harmonic  analysis,  etc. 


Maneuver  No.  of  Execution 


Case# 

Type  of  Analysis 

Time  Lgth. 
~  Seconds 

At 

Seconds 

Maneuver 

Points 

Time~ 
CP  Sec 

Total 

Cost 

Cost 
Time  Pt 

1 

TV 

TRIM+MAN.SIM. 

.201 

.00386 

52 

73.742 

$  93.82 

$1.80 

2 

QS 

TRIM+MAN.SIM. 

8.224 

.01543 

533 

455.266 

$552.60 

$1.04 

3 

QS 

TRIM+MAN  .  SIM. 

.300 

.05 

6 

27.068 

$  34.72 

$5.79 

4 

QS 

TRIM+MAN.SIM. 

2.500 

.02 

125 

125.954 

$154.28 

$1.23 

5 

QS 

TRIM+MAN.SIM. 

6.900 

.02 

345 

320.090 

$389.15 

$1.13 

16 

QS 

TRIM+MAN.SIM. 

7.500 

.02 

375 

383.842 

$479.49 

$1.28 

The  cost/time  point  for  case  3  is  ambiguous  since  it  cost  $27.89  of  the  $34.72 
total  cost  to  perform  the  trimmed  flight  analysis.  Therefore,  it  cost  $6.83  for 
6  maneuver  time  points,  $1.14/tlme  point.  It  must  be  kept  in  mind  that  in  each 
of  the  above  cases,  the  total  cost  includes  the  cost  to  perform  the  trimmed 
flight  analysis. 


From  the  data  presented  it  can  be  seen  that  it  costs  approximately 
$1.15/time  point  to  perform  a  maneuver  simulation.  Therefore,  a  maneuver  of  5  to 
10  seconds  in  length  consisting  of  250  to  500  time  points  is  going  to  cost 
between  $275  and  $575  to  run. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

1.  The  Rotorcraft  Flight  Simulation  Program  C81  provides  the  NAVAIRDEVCEN 
with  an  analytical  tool  necessary  to  evaluate  the  flying  qualities,  performance 
and  aircraft  loads  of  existing  and  future  helicopters. 

2.  The  program  in  its  present  form  here  at  NAVAIRDEVCEN  (with  the 
eigenvalue  calculation  problem  and  inability  to  perform  postprocessing 
operations)  has  a  limited  usability.  Trims  and  maneuvers  can  be  calculated; 
total  aircraft  and  rotor  partial  derivatives  can  be  calculated.  The  limitations 
are  that  characteristic  equation  roots  and  transfer  functions  have  to  be  computed 
in  another  program  and  time-history  plots  of  the  maneuver  data  must  be  hand 
plotted. 

3*  The  parameter  sweep  operation  works  as  advertised. 

4.  Due  to  the  method  of  billing  for  computer  execution  time  here  at 
NAVAIRDEVCEN,  the  cost  of  running  this  program  is  higher  than  that  which  most 
other  user  facilities  experience.  Maneuver  simulation  is  an  operational 
capability  that  is  an  important  asset  to  this  program.  It  also  is  expensive  to 
run.  It  was  shown  that  a  typical  maneuver  cost  approximately  $1. 15/time  point 
and,  at  this  rate,  it  is  not  unusual  for  a  single  maneuver  simulation  to  cost 
about  $500. 

5.  It  is  recommended  that  the  second  CDC  conversion  of  the  AGAP80  version 
of  C81  (with  postprocessing  capabilities  included)  be  installed,  checked  out  and 
validated  here  at  NAVAIRDEVCEN  as  soon  as  the  conversion  is  completed. 

6.  The  time-history  plotting  operation  is  a  highly  desirable  and  valuable 
asset  of  C81  that  can  be  used  extensively  in  conjunction  with  maneuver 
simulation. 

7.  The  eigenvalue  calculation  portion  of  the  rotorcraft  stability  analysis 
must  be  made  operational.  The  problem  is  still  under  investigation. 

8.  After  the  second  CDC  conversion  is  implemented  on  the  CDC  CYBER  175/CDC 
6600  computer  system  at  NAVAIRDEVCEN,  the  program  can  be  made  available  for 
center-wide  usage  by  other  helicopter  technology  disciplines,  i.e.,  structural 
loads  analysts,  performance  analysts,  and  control  system  analysts. 

9.  Future  Navy  helicopter  designs  can  be  modeled  and  their  flying 
qualities,  flight  loads,  and  performance  can  be  analytically  studied  with  this 
program  when  its  full  capabilities  are  brought  to  fruition  at  this  center. 


NADC-8 1290-60 


REFERENCES 


1.  Van  Gassbeek,  J.  R.,  ROTORCRAFT  FLIGHT  SIMULATION  COMPUTER  PROGRAM  C81  WITH 
DATAMAP  INTERFACE,  Volume  I  -  User's  Manual,  Bell  Helicopter  Textron, 
USAAVRADCOM-TR-80-D-38A,  Applied  Technology  Laboratory,  U.S.  Army  Research  and 
Technology  Laboratories  (AVRADCOM),  Fort  Eustls,  Virginia,  October  1981. 


2.  Etkin,  Bernard,  DYNAMICS  OF  FLIGHT,  New  York,  N.Y.,  John,  Wiley  and  Sons, 
Inc.,  1959. 


3.  McLarty,  T.  T.,  et  al. ,  ROTORCRAFT  FLIGHT  SIMULATION  WITH  COUPLED  ROTOR 
AEROELASTIC  STABILITY  ANALYSIS,  Volumes  I  -  III,  Bell  Helicopter  Textron, 

USAAMRDL  Technical  Reports  76-41A,  76-41B,  and  76-41C,  Eustls  Directorate,  U.S. 
Army  Air  Mobility  Research  and  Development  Laboratory,  Fort  Eustls,  Virginia,  May 
1977,  AD  A042464  (TR  76-41A  only). 


4.  Philbrick,  R.  D.,  and  Eubanks,  A.  L.,  OPERATIONAL  LOADS  SURVEY  -  DATA 
MANAGEMENT  SYSTEM,  Volumes  I  and  II,  USARTL  TR  78-52A,  -52B,  Applied  Technology 
Laboratory,  U.S.  Army  Research  and  Technology  Laboratories,  Fort  Eustls, 
Virginia,  1979. 


NADC-8 1290-60 


APPENDIX  A 


SAMPLE  INPUT  DATA  DECK  FOR  AH-lG  HELICOPTER 
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